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, Abstract 

The present paper is devoted to the well-posedness issue for a low-Mach number limit system with heat conduction 
but no viscosity. We will work in the framework of general Besov spaces Bp tT (M. d ), d > 2, which can be embedded 
into the class of Lipschitz functions. 

Firstly, we consider the case of p G [2, 4], with no further restrictions on the initial data. Then we tackle the 
case of any p G]l,oo], but requiring also a finite energy assumption. The extreme value p = oo can be treated 
due to a new a priori estimate for parabolic equations. At last we also briefly consider the case of any p G]l, co[ 
but with smallness condition on initial inhomogeneity. 

A continuation criterion and a lower bound for the lifespan of the solution are proved as well. In particular 
' in dimension 2, the lower bound goes to infinity as the initial density tends to a constant. 

-I— > i 

1 Introduction 

The free evolution of a compressible, effectively heat-conducting but inviscid fluid obeys the following 
equations: 

> ' 

dtp + div (pv) = 0, 

(1) { dt(pv) + div (pv <g> v) + Vp = 0, 

dt(pe) + div (pve) — div (fcVz9) +pdivv — 0, 

where p = p(t,x) G K + stands for the mass density, v — v(t,x) G W l for the velocity field and e = 
e(t, x) G K + for the internal energy per unit mass. The time variable t belongs to M + or to [0, T] and the 
space variable x is in R d with d > 2. The scalar functions p — p(t, x) and $ = i?(f , a;) denote the pressure 
and temperature respectively. The heat-conducting coefficient k = k(p,i!}) is supposed to be smooth in 
both its variables. 

| We supplement System (flj with the following two equations of states: 

(2) p = Rpd, e = Ctf, 

where R, C v denote the ideal gas constant and the specific heat capacity at constant volume, respectively. 
That is, we restrict ourselves to perfect gases. 

In this paper, we will consider highly subsonic ideal gases strictly away from vacuum, and correspond- 
ingly, we will work with the system (see (0]) or ([7]) below) which derives from System (TTJ) by letting the 
Mach number go to zero. More precisely, just as in [?], suppose (p,v,p) to be a solution of System 
(fT]) and define the dimensionless Mach number e to be the ratio of the velocity v by the reference sound 
speed. Then the rescaled triplet 

(p e (t,x) = p(^,X^, V e (t,x) = \(^,X^, p £ (t,x)=p(^,xyj 
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satisfies 



(3) 



d t pe + div(p £ v e ) = 0, 

dt(pev e ) + div ( Pe v £ ®v £ ) + ^ = 0, 

■z=i(d t pe + div (p e v e )) - div (fc e Vtf e ) + p E divv e = 0. 



Here 7 := C p /C v = 1 + R/C v represents the adiabatic index and the constant C p , denotes the specific 
heat capacity at constant pressure. The rescaled temperature and heat-conducting coefficient are given 
by 

$ £ (t, x) = , k £ (t, x) = ~k(-,Xj. 

Now let £ go to 0, that is. the pressure p £ equals to a positive constant Pq by Equations (|3]) 9 and ©3, 
and thus System ([3]) becomes formally the following zero Mach number system immediately (see [?], [?] 
and [?] for detailed computations): 

{d t p + div (pv) = 0, 

d t (pv) + div (pv <g> u) + Vn = 0, 
divv - ^div(fcVtf) = 0, 

where II = II(t , x) is a new unknown function. We remark that the justification of this limit process was 
showed in [?] in the framework of classical solutions. 

Now, again similarly as in [?], in order to make use of the almost "divergence-free" condition ((Ulg, we 
reduce System (0]) by introducing some new unknowns below. Before going into details, let us suppose 
here that the density p always has positive lower bound and converges to some constant (say "1") at 
infinity, in the sense detailed below (see (|16[) ). We also point out that in the following, the fact that 
p§ = Pq/R is a positive constant will be used thoroughly. 

Firstly for notational simplicity, we set a to be a positive constant defined by 

7-1 R 1 

a 



7P0 C p Po C pP d 

Then we define the following two coefficients, always viewed as regular functions of p: 

k = n(p) = ak{l and A = X(p) = p 1 . 

Furthermore for convenience, functions a = a(p), b = b(p) are introduced to be the antiderivatives of the 
smooth functions n(p) and —n(p)p~ l respectively, such that a(l) = 6(1) = 0. It is easy to see that 

(5) Va = KVp = aMVp = -akpVd = -pVb. 

We moreover define the new "velocity" u and the new "pressure" 7r respectively as 

(6) u = v — a/cVi? = v — V6 = v + np~ 1 \7p, tt = IT + akpdt'd = II — d t a . 
Then System (0| finally becomes 

!d t p + u ■ Vp — div (kS7 p) = 0, 
d t u +(u + V6) • Vit + AVtt = h, 
div u = 0, 

where 

(8) h(p, u) = p~ Miv (v (g) Va) = -u ■ V 2 b - (u ■ VA)Va - (V6 • VA) Va - div (V6 <g> V6) . 
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Let us just show how (JJJ 2 comes from ((4j) 2 - Observe that (J4J 2 gives 



(9) 



d t (pu) + div (pv <g> it) + dt(pVb) + div (v ® pVb) + VII = 0. 



It is easy to find that 



dt(pVb) + vn 



3 f Va + VII = V7r and div (w ® pV6) 



div (i> (g> Va). 



Thus by view of Equation (UJ) 1 , Equation © can be written as 



(10) 



pdtu + pv ■ Vu + Vir = div (i> ® Va). 



We thus multiply (dUJ) by A = p^ 1 to get Equation J7| 2 . 

Let us note that if K = 0, then a = o = and hence System ([7]) becomes the so-called density- 
dependent Euler equations 



We can see that, due to the null heat conduction, the density p satisfies a transport equation, the "source" 
term h vanishes and the velocity v itself is divergence-free. In System ((7J instead, one encounters a 
quasilinear parabolic equation for p and moreover, the transport velocity u is no longer solenoidal and 
the "source" term h is a little complicated, as it involves two derivatives of p. As early as in 1980, 
H. Beirao da Veiga and A. Valli in [?, ?, ?] have investigated the local well-posedness issue in some 
smooth bounded domain of the nonhomogeneous Euler system (fTTj) . We also cite here the book [?] as 
a good survey of the boundary-value problems in mechanics for nonhomogeneous fluids. By use of an 
energy identity, in [?] R. Danchin studied System (fTTj) in the framework of nonhomogeneous Besov space 
Bp r (M. d ) which can be embedded in C ' 1 . Recently in [?], R. Danchin and the first author treated the 
end point case where the Lebesgue exponent p in the Besov space B s p r is chosen to be 00. We notice 
that if morever p = 1 is a constant density state, then System (fTTj) reduces to the classical Euler system, 
which has been deeply studied. 

However, to our knowledge, there are not so many theoretical works devoted to the heat-conducting 
inviscid zero Mach number System @ , or equivalently System ([7]). It is interesting to view System (JJJ as 
the model for an inviscid fluid consisting of two components, both incompressible, with a diffusion effect 
obeying Fick's law. In fact, we can write Equation d3J) 3 as the following Fick's law: 



with k denoting the diffusion coefficient. Now p, u and v are considered to be the mean density, the 
mean-volume velocity and the mean-mass velocity of the mixture respectively. As usual, VII denotes 
some unknown pressure. For more physical backgrounds of this model, see [?]. One can also see [?] 
for a local existence and uniqueness result of the initial boudary value problem for this model, in the 
framework of classical solutions. In this paper, we will show some similar well-posedness results as in 
[?, ?] for the Cauchy problem for System (UJ), in the scale of Besov spaces which can be embedded into 
the class of Lipschitzian functions. 

Let us mention that the justification of the incompressible limit of Euler equations (i.e. incompressible 
Euler equations can be viewed as compressible Euler equations when the Mach number goes to zero) has 
attracted much attention. There are many early works such as [?, ?] by D.G. Ebin, [?, ?] by H. Isozaki, 
[?] by S. Klainerman and A. Majda, [?] by S. Schochet, [?] by S. Ukai, etc. for the isentropic Euler system 
case. At the beginning of this century, T. Alazard in [?] and G. Metivier and S. Schochet in [?] treated 
problems with boundary and without boundary respectivley, for the non-isentropic case. But there, the 
heat conduction did not come into play, and hence the entropy there satisfies a free transport equation. 



(11) 



dtp + v ■ Vp = 0, 
d t v + v ■ Vv + Wir = 0, 
div v = 0. 



div (v + kV In p) = 0, 
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Later T. Alazard in [?] generalized the study to various models which include the case of heat conductive 
ideal gases, i.e. he justified rigorously also the low Mach number limit from System ([T]) to System 
Let us also mention that if the fluid is viscous, that is to say there is an additional tensor term 

—diver = — div(/Lt(Vw + Dv)) — V(Vdivi)) 

in the evolution equation ^ 2 f° r the velocity v (here p, v denote viscosity coefficients which may depend 
on p), then System Q becomes the low Mach number limit system of the full Navier-Stokes system. 
See paper [?] by T. Alazard, [?] by R. Danchin and the second author, and references therein for some 
relevant results. One notices that such a viscous system can also describe viscous mixture models such 
as pollutant spreading, snow avalanche, etc.: see [?, ?, ?, ?, ?, ?, ?, ?] for some further discussion. 

Our paper is organized in the following way. 

In the next section we will present our main results on System ([7]) (the corresponding results for 
System (@| will be pointed out, in Remark 12.21 for example), i.e. the local-in-time well-posedness in Besov 
spaces Bp r c — y C 0,1 . Firstly, we focus on the case when p S [2,4] (see Theorem I2.1[) : then, we consider 
also the general instance p +00], under the additional assumptions of data in L 2 (see Theorem I2.9[) 
or when the initial density is a small perturbation of a constant state (Theorem [2T2]). We will also state 
a continuation criterion for solutions to our system, and a lower bound for the lifespan of the solution in 
terms of the norms of the initial data only: this will be done in Theorems 12.31 and 12.71 In particular in 
dimension 2, a better estimation of the lifespan under finite energy hypothesis is given in Theorem 12. Ill 

Section |3] is devoted to the tools, from Fourier analysis, which we will use in our study. In particular, 
we will briefly recall the main ideas of Littlewood-Paley decomposition and paradifferential calculus, and 
some basic properties of Besov spaces. 

In Section @] we will tackle the proof of Theorem 12.11 it will be carried out in a standard way. First 
of all, we will show a priori estimates for the linearized equations. Then, we will construct inductively a 
sequence of smooth approximated solutions. Finally, we will show its convergence to a "real solution" to 
the original equations. 

Section^ instead, is devoted to the proof of the continuation criterion (Theorem [52]) and of the lower 
bound for the lifespan of solutions to system ([7]) (Theorem I2.7[) . 

A new estimate for parabolic equations in Besov spaces r (see Proposition I6.1[) will open Section 
El This will be the key to tackle the case of finite energy initial data in this endpoint space, and so to 
prove Theorem 12.91 Moreover, in this same instance (finite energy and p = +00, s = r = 1) we will 
improve the lower bound for the lifespan of the solutions in the case of dimension 2 (see Theorem 12. lip . 
In particular, we will prove that the solutions tend to be globally defined in time whenever the initial 
density is close (in a suitable sense) to a constant state. In the same section we will also sketch the proof 
of Theorem 12.121 when the inital density is near a constant state. 

Finally, in the Appendix we will give a complete proof of the commutator estimates stated in Lemma 
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2 Main results 

Let us focus on System ([7]) to introduce our main results. In view of Equation ((TJi) of parabolic type, 
by maximum principle we can assume that the density p (if it exists on the time interval [0,T]) has the 
same positive upper and lower bounds as the initial density po : 

(12) < p* < P (t,x) < P *, Vt e [0,T], x e R d . 
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Correspondingly, the coefficients k and A can always be assumed to have positive upper and lower bounds 
too, which ensures that the pressure 7r satisfies an elliptic equation. As a matter of fact, applying operator 
"div" to Equation |Q 2 gives the following elliptic equation in divergence form: 

(13) div(AVTr) =div(/i-i>- Vu). 

Since by Equation ([7]) 2 there is no gain of regularity for the velocity u as time goes by, we suppose 
the initial divergence-free "velocity" field uq to belong to some space Bp r which can be continuously 
embedded in C 0,1 , i.e. the triplet (s,p,r) e 1 X [1, +oo] 2 only has to satisfy the following condition: 

d d 

(14) s > 1 + - , or s = l + -.r = l. 

P P 

This assumption will be enough to have the velocity field u Lipschitz continuous, and so to preserve the 
initial regularity as a transport velocity. This at least requires the "source" term h — Wtt to belong to 
L 1 ([0, T]; r ) which, by view of definition © of h and bound (|12[) for p, asks at least 

V 2 peL l {%T];B s Pir ), Vper([0,T];L°°), and Vtt g L\[0, T]; B s pr ). 

Keeping in mind that p satisfies the parabolic type equation (J7} 1 , we expect to maintain the initial 
regularity and to gain also two orders of regularity with respect to the space variable when integrating 
in time (see Subsection l4.1.1D . at least in a small time interval. We thus have to assume the difference of 
the initial density po and a constant (say 1) to be in the same space Bp r as above. However in general, 

we only get V 2 p in the time-dependent Besov space L)p(Bp r ), which is a little bit larger than L\.{Bp r ), 
see Definition 13.21 Therefore in the whole context we will deal with the spaces L^? (J3* r ) and L^^Bp r ). 

On the other hand, if h — v ■ Vu is already in space L^([0, T); Bp r ), then applying operator "A/' (see 
Section [3] for the definition) to Equation (|13l) gives 

div (AAjVvr) = div Aj(h - v • Vu) + div ([A, A^Vtt). 

So, it remains to control the commutator div ([A, Aj] Vtt), hence the low frequencies of V7T, by the classical 
commutator estimation (see Lemma 14. II for example). Due to the fact that, for Equation (|13[) above, 
we have the a priori estimate 

||Vtt|| L p < C\\h-v- Wu\\ L p 

independently of A only when p = 2 (see Lemma 2 of [?]), we have to make sure that h — v ■ Vu g L 2 . 
Hence the fact that h is composed of quadratic forms entails that p has to verify 

(15) PG[2,4]. 

To conclude, we have the following theorem, whose proof is shown in Section 01 

Theorem 2.1. Let the triplet (s,p,r) g R x [l,+oo] 2 satisfy conditions (fl"4")) and (TT5| . 

Let us take an initial density state po and an initial divergence-free velocity field uq such that 

(16) < p, < po < p* , \\po - l\\ B . ir + \\u \\ Bkr < M , 

for some positive constants p*, p* and M. 

Then there exist a positive time T (depending only on p*, p* , M,d, s,p,r) and a unique solution 
(p, u, Vtt) to System ((7]) such that (g, u, V7r) := (p — l,u, V7r) belongs to the space Ep r (T), defined 
as the set of triplet (g, u, V7r) such that 

g g C([0,T] l B; tr )nL 1 ([0,T];B p + 2 ), 

(17) I u e C([0,T];B' r ), 

t Vtt g L 1 ([0,2l;B-JnL 1 (M;£ 2 ), 
with C w ([0,T]; Bp r ) if r = +oo (see also Definition \3.2\) . 
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Remark 2.2. Let us state briefly here the corresponding wellposedness result for the original system ([4]). 
By view of the change of variables (O, we have u — Vv, Vb — Qv, where V denotes the Leray projector 
over divergence- free vector fields and Q = Id — V (see Remark 13 . 5|) . Assume Conditions (TH)) and (|15|) , 
and the initial datum (po,vo) such that 

0<P*<f>o<p*, V6(p ) - Qv Q , \\po - l|| B j r + \\Vv \\ Bkr < M. 

Then, there exist a positive time T (depending only on p*, p* , M, d, s,p, r) and a unique solution (p, v, VII) 
to System ^ such that 

p - 1 G C([Q, T);B; <r ) n Z^Q, T); B|+ 2 ) , 

« e c([o, r] ; B*- 1 ) n Z 2 ([o, r] ; s« r ) , ^ g c([o, r] ; s« r ) , 

VII G L 1 ([0,T];S* r ) , 

with C w ([0,T];^ )r ) if r = +oo. 

One notices from above that the initial velocity vq needs not to be in Bp* r but the velocity v(t) 
will be in it for almost every t G [0,T]. That is, the initial non-Lipschitzian velocity may evolve into 
a unique Lipschitzian velocity. On the other side, Theorem 12.11 doesn't say that if the initial datum 
(po — 1, vq) G r {K l ) d+1 , then there exists a unique local-in-time solution (p — 1, v, VII) to System (j4]). 
In fact, we can just say that if initially, p — 1 <E S*+ x , v G B* r , p G [p*,p*] and V6(p ) = Qvq- 

Let us also point out that since dtp L 1 ([0, T]; i 2 ) in general, we do not know whether VII G 
L 1 ([0,T];L 2 ) (recall also definition ©). Hence it seems not convenient to deal with System (Q| directly 
since the low frequences of VII can not be controlled a priori. 

Next, one can furthermore get a Beale-Kato-Majda type continuation criterion (see [?] for the original 
version) for solutions to System ([7]), similar as the ones in [?] or [?] for the density-dependent Eulcr 
System pip . As it is a coupling of a parabolic equation (for the density) and an Euler-type equation (for 
the velocity field and the pressure), one may expect, a priori, to impose conditions similar, or even weaker, 
to those found in the quoted papers, because the regularity of the density improves in time evolution. 
Actually, this is not the case: the criterion holds true under additional conditions, which are motivated 
by the structure of the nonlinearities in h and v ■ Vit. In fact, for instance, we consider L^(i?p r )-norm 
of h and V6 • Vit: 

• the term u ■ V 2 6 requires at least [|V an d IMIl~(l°°) ! 

• the term AbVb requires ||V6||/,~(£oo) and || A&Hj^ 

• the transport term Vb ■ Vu needs control on HVmH^^oo), since we only have V6 G L\(B^ r ). 
Finally, we have the following statement. 

Theorem 2.3. [Continuation Criterion] Let the triplet (s,p,r) G Kx [1,+co] 2 satisfy conditions (fl4|) 
and (|15p . Let (p, u, \7ir) be a solution of ([7]) on [0,T[ xK d such that: 

(i) p - 1 G C([0, T[ ; B s p r ) n Lj oc ([0, T[ ; B s + r 2 ) and satisfies 

(18) sup ||Vp(f)||z« + / ||V 2 p(i)||^ dt < +w; 

te[o,T[ Jo 



(ii) u G C([0, T[ ; B^ r ) and satisfies 
(19) 



sup ||u(t)||i°o + / ||Vu(t)||ioo dt < +oo; 
te[o,T[ 



G 
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(Hi) V7r G Lj oc ([0, T[ ; Bp r ) verifies, for some a > 0, 

(20) / ||V7r(t)|| fl - ;onioo dt < +00. 

Jo 

Then (p, u, VII) co-u/d be continued beyond T (if T is finite) into a solution of ([7]) with the same 
regularity. 

The proof of Theorem 12.31 issues from the following fundamental lemma whose proof can be found in 
Section 15.11 

Lemma 2.4. Lets > 0, p G (1, +oo) and r G [1,+oa]. Let (p, u, Vx) be a solution of (0 ewer [0,T[xIR d 
smc/i i/ia£ Hypothesis (i), (ii), (Hi) in Theorem \2.SH hold true. If T is finite, then one gets 

(21) \\P~ 1 k ¥ (B l , r )nL} r (B> p + 2 ) + ll U HZ§?(BJ, r ) + ll V7r (*)HiJ.(B|, r ) < +°°- 

In fact, from Theorem 12.11 one knows that once Conditions (fT4")l . (fTSf and (fTB)) hold, the initial data 
will evolve into a unique solution on a positive time interval, which depends only on the bounds on the 
density, the Bp* r -norm of the initial data and the indices. Therefore, once we have (j2"Tj) and besides 

Conditions ((Ml) . (fT5)) hold true, then there exists a time to such that, for any T < T, System ([7} with 

initial data ^p(T), u(T)\ has a unique solution until the time T + to- Thus, if we take, for instance, 

T = T — (to/2), then we get a solution until the time T+ (to/2), which is, by uniqueness, the continuation 
of (/9, u, Vtt). Theorem 12.31 follows. 

Let us analyse the scaling property of our system in order to show an explicit relationship between 
the size of the initial data and the lifespan. Recalling (|SJ), it's easy to see that, if (p, u, Vtt) is a solution 
of (|7J| with initial data (poj uq), then for any e > 

(22) (p e , u 6 , Vtt £ ) (t,x) := , e u , e 3 Vtt) (e 2 t,e x) 

is still a solution of (J7J, with initial data (pg, Uq)(x) := (po, £Mo)(£^)- So, the following result immediately 
follows. 

Proposition 2.5. Let (p, u, V7r) 6e a solution of System ([Jj wit/i initial data (po,uo) on a time interval 
[0, T[ . T/ien t/ie triplet (p e , u £ , Vir E ), given by (|22p . is sizZZ a solution of (UJ, defined at least on the 
time interval [0, T e [, wt/i 

T e > Te -2 . 

Remark 2.6. In the homogeneous case, we know that if the initial velocity is of order e then the lifespan 
is at least of order e" 1 . But here for our system (J7J (and hence for System Q), we have from the above 
proposition that if the initial data Pq — 1,Mq are of order e s ~ d /p and £ s + 1 - d /p respectively, then the 
lifespan is at least of order e~ 2 , due to the fact that 

U - l||B- r - e s - d/p \\po l||n. r , KI| B - r ~ s s - d /P +1 \\u \\ BSp r . 

In particular if s = 1 + d/p, then the lifespan is of order 0(e~ 2 ) for the initial data p§ — 1 = 0(e), 
«S = 0(e 2 ) ■ 

Even in the two-dimensional case, it's hard to expect global-in-time well-posedness for this system: 
the parabolic equation (JTJ)^ allows to improve regularity for the density term, but such a gain is (roughly 
speaking) deleted by the nonlinear term in the momentum equation (0 2 . Hence, the obstacles to global 
existence are the same as those one has to face in considering the density-dependent Euler system (|11[) . 
which has already been dealt in [?]. However, we manage to establish an explicit lower bound for the 
lifespan of the solution, depending only on the norms of the initial data, in any dimension d > 2. The 
proof will be the matter of Section 15.21 
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Theorem 2.7. Under the hypotheses of Theorem \2.1[ there exists a positive constant L, depending only 
on d, p, r, p* and p* , such that, the lifespan T of the solution to System (Q given by Theorem \ 2.1\ is 
bounded from below by the quantity 

( \ 

L 



(23) lj — , log 



J + (lIpo-ill^/IKII^J j 

where X > 2 is a constant big enough, depending only on s,d,p. 



V 



Remark 2.8. Thanks to Theorem 12. 31 the lifespan is independent of the regularity. Hence, if we want 
to get the lower bound of the lifespan, we only deal with the endpoint case s = 1 + d/p (p 6 [2,4]) and 
r = 1, whose lifespan is the largest. Therefore, the Bp r -norm in (|2"3"]) can be replaced by B^\ d ' 4 -novm. 

In Theorem 12 . 1 1 1 we will improve the previous result for 2-D fluids. Under the additional requirements 
of finite energy initial data (see also below) and taking p = +oo, we will show that the lifespan tends to 
+oo, i.e. the solution tends to be global in time, if po is "close" to (say) 1. 

As we have already pointed out, Hypothesis (fTS)) over the index p is required to get h £ L 2 , and so to 
solve the elliptic equation for the pressure term. In fact, we can remove it and consider any p S ]1, +oo] , 
provided that the following energy estimates hold: this immediately gives Vtt £ L\(L 2 ). 

First of all, if the initial datum po — 1 6 L 2 (M. d ), then by taking the L 2 (R d )-hmer product between 
p — 1 and Equation (ffil-p i.e. Equation 1(7]^. and integrating in time variable, we arrive at the following: 

i r r* r i 



(24) - / \p(t)-l\< + / / «|Vp|" - ~\\po-l\\v>. 

L JR d JO JR d z 

Such an energy identity requires the following equality to hold true: 

(u ■ Vp) (p - 1) = 0. 



In fact for instance, if u € Lj?(L°°), which is always the case under our hypothesis, then for any finite t 
it holds, by view of divu = and the previous Identity (|24l) . 

In the same way, formally, if we take L 2 -inner product between Equation (1101) and u, then thanks to 
Equation ([4"]), and integration by parts, the first two terms entail 

(25) / {pd t u + pvVu)-u = \- { P \u\ 2 . 

JRd z ore y R d 

Thanks to the divergence-free condition of u, the inner product (V7T, u) L 2^ R d^ vanishes. It rests to dealing 
with (div (v <8> Va), u) L 2^dy In fact, from Relations ([5]) and ([5]) we have 

div (v ® Va) = AbVa + u ■ V 2 a + V& • V a. 

The bound (TT2")) implies that (with C depending on p*, p*) 

HA&IU-, ||V 2 a|| L =o < C (||Vp|| 2 ^ + ||V 2 p|| L ~) and ||V6|| L 2, ||Vo||i» < C\\Vp\\ L 2. 

Thus 

(26) |(div(v® Vo),tt> za(R -)| < C-dlVpHloc + ||V 2 p|| L ~)(||Vp|| L 2 + || u || L2 )|| u || i2 

<c(e'wi| u || 2 2 + ||v P || 2 2 ), 

where we have defined 

e(*)= AllVpHi- + HVplli. + ||V 2 P |U» + |[VV||i-)dr. 
Jo 
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Hence Identity (|25|) . Holder's Inequality, Gronwall's Inequality and Energy Equality (l24|l give 

(27) |Kt)||i, < Ce c ^ (\\u \\l, + J* || Vp||i,) < Ce c ^) (|| Mo ||| 2 + || po _ i||| 2 ) . 

On the other hand, if there exists (p— 1, u) such that pT|) 1 and (flT)) ,-, hold for some finite time T > 0, 
with (s,p, r) verifying l[T3|l . then we have also (fT5|) and (TH)]| , Thus, if we assume the initial data (po, ito) 
to satisfy 

(28) Po-l,woei 2 , 

then by Energy Estimates (|24l) . (f27|). we gather it,p — 1 6 L^(L 2 ) and Vp e L^(L 2 ). Hence for any 
finite positive time T, 

ft - v ■ Vu = p^div (v <g> Va) - v ■ Vu = p" 1 Ab Va + p -1 ?; • V 2 a - u • Vu e L\(L 2 ), 

which ensures V7r S L^(L 2 ), by view of Equation (TT^J) and Lemma T3. 131 fsee below). 

To conclude, by the arguments before Theorem 12.11 and above, and a new a priori estimate for linear 
parabolic equations in Besov spaces B^ r (see Proposition 16. II) . we have the following statement. 

Theorem 2.9. Let the triplet (s,p,r) € M x [l,+oo] 2 satisfy Condition (TUT) and 1 < p < oo. 

Let the initial density po and the divergence-free initial velocity uq fulfill (|28p as well as (|16p . /or some 

positive constants p*, p* and M . 

Then there exist a positive time T (depending only on p», p* , M, d, s,p,r) and a unique solution 
(p,U, Vtt) to System © with (p-l,u, Vvr) e E» r (T) such that p-l,u£ C([0,T];L 2 ) andVp e L^(L 2 ). 

Remark 2.10. Notice that we have the embedding £?^ r (K d ) i 2 (R d ) if s > d/p - d/2, p 6 [1, 2] (see 
Corollary [X71 below) . Thus if the triplet (s,p,r) e K x [1,+oc] 2 satisfy Condition (fTI]) and 1 < p < 2, 
then the initial data po,uo satisfying (TTB1) alreday verifies finite-energy condition (j2"5)) . 

In fact, if we take p = ~p constant in system ([7J, or equivalently in (jj), we get the classical homogeneous 
Euler system. For this system, the global-in-time existence issue in dimension d — 2 has been well-known 
since 1933, due to the pioneering work [?] by Wolibner. For non-homogeneous perfect fluids, see system 
(fTTj) , it's still open if its solutions exist globally in time. However, in [?] it's proved that, for initial densities 
close to a constant state, the lifespan of the corresponding solutions tends to +oo (independently of the 
initial velocity field). 

Actually, under the same hypothesis of Theorem 12.91 in dimension 2, we are able to prove a similar 
result also for our system. 

Theorem 2.11. Let d = 2. 

Let us suppose the initial data po and uq to be such that go := po — 1 G L 2 1 , with < < po < p* , 

and u 6 L 2 H £^,1- 

Then the lifespan of the solution to system ([7]), given by Theorem \2.9l is bounded from below by the 
quantity 

<29) l + IMIi. + Vll^., bS (' + bS ((l + IW 2L)Mlej ) ' 

for a suitable exponent DJl > 5 and a constant L which depends only on p* and p* . 

In particular, the lifespan tends to +oo, i.e. the solution tends to be global, for initial densities close 
(in the B^ 1 norm) to a constant state. 

There is another way to get rid of Condition (|15[) imposed on p, that is, we assume an additional 
smallness condition over the initial density, which ensures that the elliptic equation for the pressure is 
almost the Laplace equation, up to a perturbation term. 
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Theorem 2.12. A small constant c > exists such that the following statement holds true. 

Let the triplet (s,p,r) £ K x [l,+oo] 2 satisfy Condition (|14p and p 6 (l,+oo). Le£ ws take an initial 
density state po and an initial divergence-free velocity field uq such that the bounds in (|16p are true, for 
some positive constants p*, p* and M. Assume moreover \\po — < c. 

Then there exist a positive time T (depending on p*, p* , c, M, d, s,p,r) and a unique solution (g, u, V7r) 
to System belonging to the space I£ r (T) , defined as Ep* r (T) ( recall (|17p ), but without the condition 
Vtt G i 1 ([0, T]; L 2 ). 

Remark 2.13. It goes without saying that under hypotheses of Theorem 12.91 or Theorem 12. 121 a con- 
tinuation criterion and a lifespan lower bound analogous to those of Theorem 12.31 and Theorem 12.71 
respectively, can be proved. The corresponding results for the original system (fj[]) also hold true, similar 
as in Remark 12.21 which are omitted here. 

We change the word here that in the sequel, C always denotes some "harmless" constant which may 
only depend on d, s,p,r, p», p* , unless otherwise defined. Notation A < B means A < CB and A ~ B 
says A equals to B, up to a constant factor. For notational convenience, the notation g always represent 
p — 1, unless otherwise specified. 

3 An overview on Fourier analysis techniques 

Our results mostly rely on Fourier analysis methods which are based on a nonhomogeneous dyadic 
partition of unity with respect to Fourier variable, the so-called Littlewood-Paley decomposition. Unless 
otherwise specified, all the results which are presented in this section are proved in [?]. 

In order to define a Littlewood-Paley decomposition, fix a smooth radial function \ supported in (say) 
the ball 5(0, |), equals to 1 in a neighborhood of 5(0, |) and such that r H> x( r e r) is nonincreasing over 
R+, and set <p(£) = *(§) - X (0- 

The dyadic blocks (Aj)jgz are defined bjQ 

Aj := if j < -2, A_i := x{D) and A,- := ip{2~ j D) if j > 0. 
We also introduce the following low frequency cut-off: 

SjU := x(2~ j D) = A i' for 3 - °> S i u = for 3 - °' 

The following classical properties will be used freely throughout the paper: 

• for any u S S 1 , the equality u = ■ AjU holds true in S'\ 

• for all u and v in 5', the sequence (Sj—iu Aju)jgN is spectrally supported in dyadic annuli. 
One can now define what a Besov space B^ r is: 

Definition 3.1. Let u be a tempered distribution, s a real number, and 1 < p, r < oo. We set 

i 

h\\B lr ~(^2 r J s \\A jU \\ r Lp y ifr<oo and \\u\\ B ^ := s^ip (2^ || A jU \\ LP ) . 

We then define the space Bp r as the subset of distributions ue5' such that r is finite. 

When solving evolutionary PDEs, it is natural to use spaces of type Lj,(X) — L p (0,T; X) with X 
denoting some Banach space. In our case, X will be a Besov space so that we will have to localize the 
equations by Littlewood-Paley decomposition. This will provide us with estimates of the Lebesgue norm 
of each dyadic block before performing integration in time. This leads to the following definition: 

throughout we agree that f(D) stands for the pseudo-differential operator u h- > T~ 1 (/(f)^ r ' u (£))- 
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Definition 3.2. For set, (q,p,r) G [l,+oo] 3 and T G [0, +00], we set 

\A jU (t)\\ q LP dt) 




with the usual change if r = +00 or q = +00. 

We also set C T {B s pr ) = r ) D C([0,T]; fl* p ). 

Remark 3.3. From the above definition, it is easy to show that for all s £ R, the Besov space Z?| 2 
coincides with the nonhomogeneous Sobolev space H s . Let us also point out that for any k G N and 
p G [1, 00], we have the following chain of continuous embedding: 

where iy fe ' p denotes the set of L p functions with derivatives up to order k in L p . The following embedding 
is also true: 

L q T { B lr) ^ L q T {B s p , r ) if 9 > r, 2|(^ >r ) l£(B* r ) if g < r . 

The Besov spaces have many nice properties which will be recalled throughout the paper whenever 
they are needed. For the time being, let us just recall that if Condition (|14p holds true then B p r is an 
algebra continuously embedded in the set C 0,1 of bounded Lipschitz functions (see e.g. [?], Chap. 2), 
and that the gradient operator maps B p r in B^" 1 . The following result will be also needed: 

Proposition 3.4. Let F : M. d — > R be a smooth homogeneous function of degree m away from a neigh- 
borhood of the origin. Then for all (p, r) G [1, 00] 2 and s G R, the operator F{D) maps B't r in B 



s—m 



Remark 3.5. Let V be the Leray projector over divergence free vector fields and Q := Id — V . Recall 
that in Fourier variables, we have for all vector field u 

Vu(0 = J2(5 jk + 1) |p u k (0 and Qu(0 = - X £ • «(£). 

Therefore, both (Id — A-i)V and (Id — A_i)Q satisfy the assumptions of the above proposition with 
m = hence are self- maps on B p r for any s G R and 1 < p, r < 00. 

The following lemma (referred in what follows as Bernstein's inequalities) describes the way derivatives 
act on spectrally localized functions. 

Lemma 3.6. Let < r < R. A constant C exists so that, for any nonnegative integer k, any couple (p, q) 
in [1, 00] 2 with p < q and any function u of LP , we have for all A > 0, 

SuppuGB(0,Ai?) => \\V k u\\ Lq <C k+1 X k+d( ^^\\u\\ LP ; 
Suppw G {£ G R N I r\ < |e| < R\} => C- k - 1 X k \\u\\ LP < \\V k u\\ LP < C k+1 X k \\u\\ LP . 

The first Bernstein inequality entails the following embedding result, which is a generalization of 
Remark 13.31 



Corollary 3.7. The space B^ ri is continuously embedded in the space r2 whenever 1 < p\ < P2 < 00 
and 

S2 < s\ — d/pi + d/p2 or S2 = S\ — d/p\ + d/p2 and 1 < r\ < ^2 < 00. 

Let us now recall the so-called Bony's decomposition introduced in [?] for the products. Formally, 
any product of two tempered distributions u and v, may be decomposed into 

(30) uv = T u v + T v u + R(u, v) 
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with 

T u v := y~^,Sj_ iuAjV and R(u,v) :='S~^ AjuAj>v. 

3 3 |i'-j|<l 

The above operator T is called "paraproduct" operator whereas R is called "remainder" operator. 

We have the following classical estimates in Besov spaces for the paraproduct and remainder operators. 

Proposition 3.8. Let s, si, S2 6i, 1 < r, r\, r2,p < oo with i < min{l, — + — }■ 

• For the paraproduct we have the following two separate estimates: 

\\T u v\\ bS1+S2 ^ <\\u\\ B;i Jv\\ B;2r2 if Sl <f, 
P '\\T u v\\ B s r <\\u\M\v\\ Blr - 

• For the remainder, if S\ + S2 + dmin{0, 1 — -} > 0, then 

\\R(u,v)\\ gSi+S2 ^ <\\u\\ B; iJv\\ B; 2 r2 . 

By using (|30[) . but keeping the operation pertaining to the time variable before taking the ^ r -norm, 
we get the following estimate concerning time-dependent Besov spaces, which will be of constant use in 
this paper. 

Proposition 3.9. For any (s,p,r,q,qi,q2) £ R x [l,oo] 4 with 

. ,d d, d , d 111 

— mm{ — , — } < s < —. with r = 1 if s — — , and — = 1 , 

VP V P q 9i 92 

£/ie following holds: 

(31) IIHIz^s^.) ^ c, ll«llzft(B^»)H' , llz«(fl; if .)- 

Moreover, if s > 0, /or any g,^ 6 [l,oo] mi/i i = 1,...,4 smc/i t/ia£ 1/g := 1/qi + X/q% = 1/93 + 1/94; 
the following also holds: 

IIHIz«(B| r ) < ^(iNz^Ci^lMlz^CBj,.) + ll"lli*(B; ir )ll' , llL«f(i-))- 

The following results pertain to the composition of functions in Besov spaces: they will be needed 
for estimating functions depending on the density and coming into play in our equations. The proof of 
Proposition 13.101 can be found in [?], while for Proposition 13 . 1 II we refer to [?]. 

Proposition 3.10. Let f be a smooth function such that /(0) = 0, and take s > and {q,p, r) £ [1, +oo] 3 . 
For any u £ Cq° , the following inequalities hold true: 

\\fou\\ B}r < C\\u\\ B s r , \\f°u\\l* T(B s r) < C\\u\\ tUB; ^, 
where the constant C depends on s,p,r, f and \\u\\ z,<»(z,oo). 

Proposition 3.11. Let I be an open interval of R and F : I — > R a smooth function. Then for all 
compact subset J C /, s > and (g,p, r) £ [l,+oo] 3 there exists a constant C such that for all function 
a valued in J and with gradient in L^(Bp^ r 1 ) 1 we have V(-F(a)) £ Lj,(B^~ ) and 

ll v ™)llz^( S -)<Gllv a |lz^ (srrl) . 
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In the analysis of ([7]) 2 , we will need a priori estimates in Besov spaces for the transport equation 

(T) 



dta + v ■ Va = /, 
a\t=o — a a- 



These estimates are provided by the following proposition, whose proof may be found in e.g. [?] 
Proposition 3.12. Let (p,r) G [l,+oo] 2 and 



(d d) . (d d) .„ ,. 

< — , — > or c > — 1 — mm < — , — > it div v = . 
{p p> } {p p' } 



(32) a > -mini-,— 
with strict inequality if r < +oo. 

Let do G Bp r , f G L 1 ([0, T]; i?p r ) and v be a time dependent vector field in L q ([0, T}; -B^/oo) f or some 
q > 1 and M > such that 

Vv G L 1 ([0,r];Bp, / t £ nL°°) if a<l + d/p, 

Vv G L^CT]^;; 1 ) if a>l + d/p, or <r = 1 + d/p and r = 1. 

T/ien Equation (T) /ias a unique solution a in 
• the space C([0,T]; -Bp r ) if r < oo, 

. fte .space (fV <CT C([0, T]; B^j) n^CMj^oo) if r = oo. 
Moreover, for all t G [0,T], we have 

(33) e- cv -W||a(t)|| Sp% < ||a || Bp% + /' e~ Cv ^ \\f(t') \\ B; p dt' 



wii/i V"(t) := 



|V«(t)||^/ £nLDO ?/ cr < 1 + d/p, 

|Vw(t)|| BCT -i if o~>l + d/p, or a = 1 + d/p and r = 1, 



and, if equality holds in (13211 and r = oo, V'(t) = ||Vi>|| n «t/j>. 

If a = v then, for all a > (cr > —1 i/divv = 0/, Estimate (I33p ao/ds imt/i V'(t) := || Va(i) H^oo . 
Finally, we shall make an extensive use of energy estimates for the following elliptic equation: 

(34) - div (a VII) = div F in R d 
where a — a(x) is a given suitably smooth bounded function satisfying 

(35) a„ := inf a(x) > 0. 

We shall use the following result (see the proof in e.g. [?]): 

Lemma 3.13. For all vector field F with coefficients in L 2 , there exists a tempered distribution II, unique 
up to constant functions, such that VII G L 2 and Equation (|34|) is satisfied. In addition, we have 

(36) a*||VII|| L2 < \\F\\ L 2. 
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4 Proof of Theorem 12.1 



The aim of this section is proving our first well-posedness result, i.e. Theorem 12. 1[ for System 0. The 
a priori estimates we establish here will be the basis throughtout the following context. 

The notation fj will always denote Ajf, unless otherwise specified. We also notice here that by 
embedding results stated in Remark l3.3l and Corollarv l3.7l for any e > 0, the following chain of embeddings 
holds true: 



which will be frequently used in our computations. 

4.1 Linearized equations 

In this subsection we want to establish a priori estimates for the linearized equations associated to original 
System ([7]). Subsection 14.1.11 deals with the linearized equation for the density. We first need two lemmas 
(see Lemma FTTl and Lemma |4~2"1 below) . corresponding to two types of commutator estimates respectively. 
Subsection 14.1.21 is dedicated to the velocity equation, where we shall pay attention to the pressure term, 
whose low frequencies can be controlled by a direct use of Holder's Inequality. 

Let us point out here that in order to prove the uniqueness, we will consider the difference of two 
solutions to System ([7]). Noticing that Equation ([7]) 2 f° r the velocity u will cause one derivative loss, we 
therefore also have to look for a priori estimates for the unknown in Bp r , under a weaker condition on 
the indices (see (|4"Tj) below), rather than flT4"l) . One notices also in this case, the "given" transport velocity 
field w in the linear equation for u (see (1551) below) is of higher regularity than u, and hence we will use 
throughly the following equality (which holds due to div it = 0): 

(38) div (w ■ Vu) = div (u ■ Vu> — u div w) . 

In addition, in proving a continuation criterion, we only need L°°-norm instead of -norm (with 
(s, r) satisfying (|14p ). Therefore, in the present paragraph we will give also a priori estimates involving 
L°° norm, i.e. we will not only make a rough use of the embedding B^ 1 <—t L°°. 

4.1.1 A priori estimates for the density 

First of all, we consider the linearized equation for the density term, which is actually the same for both 
p and q. Let us suppose the right hand side of the first equation of ((7) to be some scalar function /, a 
more general case which turns out to be useful in the sequel: 



This is a parabolic equation and the treatement of the transport term u • Vp is very classical. In fact, 
after applying the operator Aj, it rests to handle the commutator [u, A.j] ■ Vg, for which we have the 
following lemma (see [?], Chapter 2). 

Lemma 4.1. Let us take the triple (s,p,r) verifying 



(37) 



(39) 



dtQ + u ■ — div (k Vq) = f . 



(40) 



s > — c?min{-, — } and (p, r) G [1, +00] 2 , 



with 



r = 1 if s = 1 H — . 

P 



Then we have (with some positive constant C depending only on d,s,p,r) 



(41) 



I 




1 dr, 
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where {& s } s *s defined by 
(42) $ s (t) = 



l|v*>(t)|| B . 

Moreover, if s > we can also infer 



B p,i 



otherwise. 



(43) 



V s \\[^A 3 }^\\ LP \\ er dr < C f (||V^|U ||^|| B . r + HV^II*.-! HVVIIj 



dr. 



On the other side, the term div (kVq) gives us the commutator div ([k, Aj]Vg). We can do similarly 
as in [?]. But since we work in the Besov space r with r not necessarily being 1, we have to resort to 
a new commutator estimate to get the norm || • Hjiq (instead of || • ||lJ( )) °f Qi which will be absorbed 
by the left hand side. Compared with Lemma I4TT1 this is nothing but taking the interpolation before the 
integral with respect to the time variable. More precisely, we have the following lemma, whose proof can 
be found in Appendix. 

Lemma 4.2. Let the triple (a + l,p, r) satisfy (|40| with e > and a\ < a + 1 < o~2- Define 8 € ]0, 1] 
such that 

(X+l = (9(7! + (l-6») cr 2 . 

Then, for any space derivative dk, with 1 < k < d, we have (with some positive constant C depending 
only on d, o~,p,r) 



(44) 



V° I \\d k ([^A 3 }Vij)\\ LP dr 

ce 



< 



with defined by 

(45) $ a+1 (t) 



l|Vp(f)|| i i 



L f ( B p,.) 



if in addition Ui, i = 1,2, satisfy 



d d 
Oi < 1 + - if a + l < 1 + - 
P P 



z/ae(-l-dmin{i,i},|), 



if a = - with r = 1 or a > -, 



d d 

(Xi > 1 + - if (7 + 1 > 1 + - 



// moreover a > and, /or some £1, (2 fl^d 77 € ]0, 1] , we raave also a + l = rj £1 + (1 — rj) £2, i/ien 



2* / |j5 fc ([^,A,]V^)|| iP dr 



< 



(46) 



Next lemma is in the same spirit of the previous one, by view of Proposition 13 . 91 which gives estimates 
for the product of two functions (see again the Appendix for the proof). 

Lemma 4.3. Let a > and (p,r) £ [l,+oo] 2 . Fix e > and cri < a < 02. Define 8 g]0, 1] such that 
a = 8a 1 + {1-8)0-2. 
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There exists a constant C , depending only on d, a , p and r, such that, for all functions /, g G <S, one 

has 

Wf9kl(B« r) < ^y-e [ (ll/ltf- H^IIb- + \\9\\)H II/HbR.) *r 

+ (l-0)e (||/||£x (flp?r) + \\9\\lUb;%)) ■ 

Remark 4.4. Let us point out that the previous results can be also generalized in the following sense. 
For instance, as in commutator's estimate (|4"6")h take this time 9 and r\ G ]0, 1] such that 

a = 9 a x + (l-0)o-2 = vCi + (1 — »|) Ca , 
with again cri < a% and £i < C2- Then we have, for any e and 5 € ]0, 1[ , 

w *■* Q 

+ (l-f)e|l/[[£j(^? r ) + ^-^dhki^y 

for a constant C which depends only on the dimension d and on the indices a, p and r. 

Let us come back to the linear equation . Thanks to inequality (l4"Tj) and (|4"6"]l in Lemmas 14.11 and 
14.21 the following a priori estimates hold true, similar to Proposition 4.1 in [?]. 

Proposition 4.5. Let the triple (s,p,r) verify 

(47) s > — , p Gil, +oo[ , r € [1, +oo] , wit/i r = 1 if s = 1 H — or — . 

/' /' /' 

Lei u be a smooth divergence-free vector field such that Vu € B^ 1 , and k a smooth real function such 

that k > ft* > ; Vft € Moreover, let us assume the external force f € L^, (.B*). 

ITien i/iere exists a positive constant C\ depending only on n*,d, s,p,r such that, for every smooth 
solution g of (|39p with initial condition g\ t =o = Qa, the following estimate holds true for every t € [0,Tq\: 

( 48 ) \\ehf(B. ir ) + WekuBZ?) < C,e C ^ (\\Q Q \\B lr + II/IIl^)) . 
where we have defined K(0) = and 

(49) K'(t) := 1 + ||V«|| i + ||V«|| B . r + llVKlli-. 

B^nB}^ 

Proof. We follow the proof of Proposition 4.1 in [?]: we apply the operator Aj to the equation, we integrate 
first in space and then in time; then we use the commutator estimates and Gronwall's Inequality to get 
the result. 

Applying the localization operator Aj to Equation (p?9"|) yields 

(50) d t Qj + u ■ X7q 3 - div {kVqj) = fj + U) - K) , 

where we have set Qj :— Ajq, fj := Ajf, TZh :— [u, Aj] ■ V g and 1Z? := div [ft, Aj} Vg. 

According to the following Bernstein type inequality given in Appendix B of [?] (which holds true for 
any p e (1, oo) and j > 0) 

(51) -/ diy( K Vg j )\g j f- 2 Qj = (p-l) [ «|V&|%r 2 

JR d JR d 

>C(d,p,K*)2 2 > [ \g,\P, 
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and taking the L 2 scalar product between (JSOJ) and \gj\ p 2 Qj gives us, for all j > 0, 

d '\ ej \\l P + C2 2 > \\ Qj \\l P < cwe.wi- 1 (\\fj\\ LP + \\n]\\ LP + \\n 2 \\ LP ) . 



dt 

Now, we divide both members of the previous inequality by HgjH^p 1 and we integrate in time; then 
we multiply by 2 JS and take the l r norm with respect to j. Hence, for all t G [0, To] we get: 



l£'IIZ t 1 (B=+ 2 



< C[ [[fio||fl. r + 2-(^+ 2 ) \\±- X Q\\ LULf) + 



(52) 



2JS /VlWlL^ 7 " 
Jo 



+ 


V s [ 




Jo 



The low-frequency term A_ig can be easily bounded in [0,To]. As a matter of fact, by definition we 
immediately have, for all t, 



(53) 2-(*+ 2 ) \\A. ie \\ Ll{LP) < C JjQh 

For the first commutator term, by Lemma 14. 11 we get 



dr. 



(54) 



2js \\n\ LP 



dr < 



|V«|| a \\g\\ 



p.i 



ifse(-dmin{§,±},l + §), 



|Vu|| B -i Hells- p if * > 1 + i , and r = 1 if s = 1 + &. 



For the second commutator term, instead, we apply (1461) with a = ax = s, (7% = s + 2 and 77 = 1 , 
£1 = s + 1, and we get, for any small s > and some corresponding constant C e , 



(55) 



2 JS 



/ ML* dr £ £ IMb(B5t 2 ) + C 's f (llV^lli^ + ||V K || 2 Loo ) ||e|| B . r dr 
Jo Jo v ' 



for all s > d/p, with r = 1 if s = d/p, such that embedding Bpf^ L°° holds true. 

We put (|53p . and ([53)1 into |5^)1 . choose sufficiently small e and apply Gronwall's lemma to arrive 
at (gU). □ 



4.1.2 The linearized equation for velocity field and pressure term 

The linearized equation for the velocity reads 



(56) 



dtu + w ■ Vit + AV7r = ft, 
div u = 0, 
u|t=o = Uo, 



where the initial datum uq, the transport vector field w, the coefficient A and the source term h are all 
smooth and decrease rapidly at infinity. We have the following result. 



Proposition 4.6. Let 

(57) s>--~, pe[2,4], re[l,+oo], 

p 4 



d d 
lOTfft r = 1 ij s = — or H — . 



Suppose < A* < A(i, x) for all t > and x G K d . 
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Then there exists a positive constant C2 such that for any smooth solution u of (|56p . the following 
estimates hold true: 

(58) IMIZr(B |ir ) < ^^W(||^|| fl . f . + ||/ l [[ Ztl(B . (r)nii(ia) ) ) 

(59) HV7r|| £ i (fl5ir)ni i (La) < C2(||^||£i (B ^ )nL j (L2) + W(t)||M|| Zr( ^ r) ), 
where W{§) — and 



(60) W'{t) = 



\\vw{t)\\ „ ifs e(--^i + -}, 

B* a P 4 p 

||Vw(t)|| s ,-i ifs>l + -, 

p ' r p 



and the constant Ci depends on d,p, s,r, A*, A* with 



(61) \*(t) := HAlU-^oo) + 



IIVAL * i fse(--j,l + -} 



Proof. From Proposition 13 . 1 21 we easily get the following estimate for u with W(t) defined by (jSP")) : 

(62) K*)Hrr(B 5ir ) ^ eCfes)M/(t) (ll w ollB^ + ||A - AV7r|| £i(B . r) ). 

By decomposing AV7r into 



(A_iA)Vtt + [{Id - A_i)AJ Vtt, 
and by product estimates given in Proposition 13.91 we have 

(63) ||AV7r|| £tl(B » r) <CA*||V7r||j KB » r) , if s > - mm{-, 4>, with r = 1 if s = -, 

with A* defined by (foTj) . Hence it is sufficient to estimate ||V7r||ji^ BS y in order to obtain Estimate lj58|). 

By view of Equality (f3"8")h applying the operator "div" to the first equation of ([55)) yields the elliptic 
equation for 7r: 

(64) div (AV7r) = div (h — w ■ Vu) = div {h — u ■ Vw + u div u>) . 

Similar as to get (|S"2")) . we apply the localization operator Aj to we multiply \nj\ p ^ 2 TTj 1 then we 

intergrate with respect to space variable and we use the Bernstein type Inequality ([5"Tjl and Holder's 
Inequality and we finally find 

A* 2^11^11^ < HAj-div (h- w -Vit) \\ L p + ||div [A, Aj-]V7r||i», Vj > 0. 

On the other hand, for p > 2, Bernstein's Inequalities given in Lemma 13.61 entails 

||A_ 1 V7r|| i p < C||W|| L2 , 2^11^-11 ~ ||V7r,-||, Vj > 0, 

which ensure that 

(65) ^*l|V7r|| £ i (B , j < A*||V7r|| it i (i2) + \\div (h - w • VuJUji^.-ij 

+ ||2-''( s - 1 )||div[A,A,]V 7 r|| ifl(iP) ||^. 
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Let us focus on the above commutator term for a while. Commutator Estimate (|44[) with 6 = 1 in 
Lemma T4.2I implies that it can be controlled by 

||^- 1 )||div[A ) A J ]V7r|Uj (2>) || < r < A*||V7r|| £i(B .- rl)) if a > -mm{^}. 

Motivatived by the embeddings (1571) and 



L 2 S p % 2 , Vp > 2, 



we apply interpolation inequalities between Besov spaces, obtaining 

for some e s (0, 1), such that 

s + e - 1 - (d/p - d/2) , . d d 

" = s-Wp '-d/2) €(Q ' 1) ' ^ S> ^2 +1 - £ - 
Therefore, one can use directly Young's Inequality on the above bound for commutator term 
||div [A, A,] Vn\\ Ll(LP) \\ t r < C(d, s,p, r, e, e)(A*) V(i-fl) || V7r|| L i (ia) + £ || VttH^^, 
such that for some small enough e > 0, Estimate (|65|) becomes 

(66) II V7r|| zj < C7(d, a.p, r, e, A.XA*) 1 /^-^ || V7r|U i(i >) + CHdiv (/» - to - V«) Ilzj^-X). 
Thanks to Lemma T3. 131 we have already got, by Equation (l64l) . 

HVttHji^^ < C(d,s,p,r,\*,\*)\\h-u- Vw + udivw|| L i (L2) 
+ C||div (h-wVu)\\ iUB s- ly 

It rests to dealing with 

II"" • Vw|| L i (i 2), ||udivw|| L i (L 2) and ||div (w ■ Vu)\\ Z i( B s-i y 
We can easily find, by Holder's Inequality and embedding results, that 

ll M ■ v HIl1(l 2 ) < J IMM|Vtt;||i4 < J II^IIb*^ II V^H^^ dr, Vsi,s 2 >0. 

Hence for p < 4, s > — — ^ we have 
p 4 

(67) ||«-V«|| i}(i2 )< f ||«|| fl;ir ||V«>[| & dr. 

The term wdiv w is actually analogous to the previous one. 

On the other hand, recalling the divergence-free condition over it, it is easy to decompose ||div(w 
Vu)|| B s-i into 

(68) WT^djv* + T djU .d lW i + div (R{w s ,d j u))\\ B .-i, 
which can be controlled, according to Proposition ^. 8[ by 



W'(t)\\Vu\\ B ,-i, if s >-min{-,4}, 

with W'(t) defined by (JSQJ. 

To conclude, in the case j) g [2,4], s > j - |, estimate ([59")) holds and so does estimation (1551) . by 
view of (|62p and Gronwall's Inequality. 

□ 
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4.2 Proof of the existence 

In this section we will follow the standard procedure to prove the local existence of the solution to 
System ([7]): we construct a sequence of approximate solutions, we show uniform bounds and we prove 
the convergence to a unique solution. We note here some related key points: 

• Since we admit also large initial density po, we will introduce the large linear part pz, of the solution 
p as the solution to the free heat equation with the same initial datum, which is explicit and has 
positive lower bound. The remainder part p := p — pi, is small and hence easier to handle. 

• In the convergence part, since there are quantities like Vm™ € appearing in the source term of 
the equation for the "difference" sequence 8u n , we first show that the built sequence converges to a 
solution in a space with lower regularity (i.e. in space Ep^r(T), see (TTTJ)) and hence in the desired 
Besov space because of the uniform bounds for the sequence. 

• Whenever the indices s,r satisfy s>-ors>-.r — 1. the inequalities 
(69) 

\\uv\\ B}<r < C(d,s, P ,r)\\u\\ B} Jv\\ B}tT , \\f(p)\\ B . >r < C(f, IMMUfjII^, with f(l) = 0, 
and their time-dependent version 

\\uv\\ I?(B , ) ^ c {d,s,p,r)\\u\\z< 1(B , Jwllz^fB, ), with i = — + — , 

ll/(p)llzf( B . >p) < c(f, ||p|U»(i»))||e|li. (B . r j, with /(i) = o, 

always hold true, thanks to Propositions 13.91 and 13.101 

In the coming proof, C always denotes some "harmless" constant which may depend only on p*, p* , d, s,p, r 
and for simplicity we assume = A(l) = 1. 



4.2.1 Step 1 Construction of a sequence of approximate solutions 

In this step, we take (s,p,r) such that Conditions (TT4")) and (TT5|) hold true. Let us introduce the approx- 
imate solution sequence {(g n , u n , V7r ra )} n >o by induction. 
Without loss of generality we can assume 

(70) ^ < S nPo , Vn G N ; 

then, first of all we set (g° , u° , \7tt°) := (Sogo, SqUq, 0). Let us note that these functions are smooth and 
fast decaying at infinity. 

Now, we assume by induction that the triplet (f» n— u™ -1 , V7r n_1 ) of smooth and fast decaying func- 
tions has been constructed. Let us suppose also that there exists a sufficiently small parameter r (to be 
determined later) , a positive time T* (which may depend on r) and a positive constant Cm (which may 
depend on M) such that 

(71) ^<P n - 1 --=l + 9 n -\ ll^-lr-^^M, \\Q n -% l , iBi p )n Z^ B ^ ) <r, 

(72) K-ij^ < C M , IK-i^.^,^) < r, IIV^- 1 !!^^,^^ < . 

Let us immediately remark that the above estimates (|7Tj) and (|T2"|) obviously hold true for (g°, u°, V7r°), 
if T* is assumed small enough. 

Now we define (g n , u n , V7r n ) as the unique smooth global solution of the linear system 



(73) 



d t g n + u n - x ■ \7g n - div (^^Vg") = 0, 
d t u n + (u"- 1 + V6™- 1 ) • Vu™ + A n - 1 V7r n = h^ 1 , 
divu" = 0, 
{g n ,u n )\ t= o = (S n g ,S n u ), 
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where we have set 

b n ~ x = bip 11 - 1 ) , k 11 - 1 = k^- 1 ) , A"" 1 = X(p n ~ 1 ) , h 11 - 1 = ft^"- 1 ,^- 1 ) . 
We want to show that also the triplet [g n ,u n , Vtt") verifies (J7TJ) and ([72]). 

First of all, keeping in mind f|70[) , we apply the maximum principle to the linear parabolic equation 
for g n , yielding p n := 1 + g n > p*/2. 

Now, we introduce as the solution of the heat equation with the initial datum go £ Bp r : 

( d t gh - Ag L = 
\ (ei)| t =o = So- 

Then, it's easy to see (e.g. applying Proposition I4.5[) that the global solution g^ = e tA go satisfies, for 
any positive time T < +oo and some constant Ct > depending on T, 

(74) \\Q L \\z ¥i B} <r ) + WeLh^BX?) ^ CtUoWb;,,- 

We claim that, given r > 0, we can choose T* < +oo such that one has 

( 75 ) ll^llz?,»(Bj+ 1 )niJ,»(B|+ 2 ) ^ r2 - 
Indeed, we can write 

rT* \ 



QL 7. 



2i(*+ 2 ) 



e tA A,-eo|| L pdt 



Now, we use the fact that (see e.g. Chapter 2 of [?] for the proof) for any j > —1, the operator e tA Aj 
belongs to C(L P ) :— {A : L p L p linear and bounded}. Furthermore, for j > 0, its norm can be 
bounded in the following way: 

II tA a || < r -Ct2 2 i 

On the other side, we can decompose go into low-frequency, large part go,z and high-frequency, small part 
go, hi such that, for some N large enough, 

Qo,l = do on 2 N B, g j = outside 2 N B, g 0> i + g ,h = go, \\Qo,h\\B* r ^ r3 ' 
where B is the unitary ball centered at the origin. Therefore, once one notices that 



T* 



2^e- ct22j dt 



< C(l - e 



), 



2^e- ct22j dt 



j>N 



we get 



< C(l - 



-C2 2 "T' 



-C2 2 "T* 



+ c 



(2^||A^ || iP ) 



j>N 



\Qo\ 



So one can choose sufficiently small T* such that (^ s + 2 ) — r2 - The term ||gi||j2 can ^ e 

handled in the same way or by interpolation inequality. Hence, our claim (|75|) is proved. 

Now we define the sequence g^ — S n gL'- it too solves the free heat equation, but with initial data 
S n go. Hence, it too satisfies ([74")) and (|75)) . for some T* > independent of n. 

We next consider the small remainder g n := g n - g n h . We claim that it fulfills, for all n e N, 



(76) 
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d t g n + u 71 - 1 ■ X7g n - div (n^Vg 71 ) = -u"" 1 • V Q n L + div ((k™" 1 - l)Vp£), 

r|i=o = 0. 



In fact, g n = Q n — qH solves 
(77) 

So, if we define 

K n ~ l {t) := t+ || Vu™- 1 
by (gSJ) we infer that 



< Ce cr " 1(r) |-« n - 1 .Ve2 + div((/s n - 1 -i)Ve2 



Now, inductive assumptions ([71]) and (172"]) tell us that if" X (T*) < C r if T* < r and r is small enough. 
One also refers to Proposition 13.91 and Proposition 13.101 to bound 



u n - L - VqI + div [(K n - L -l)Vg 



by the quantity (up to a constant factor) 

and hence by Cm t 2 . Therefore, (1751) is proved, and hence (fTTj) holds for g" = + g£ i f° r sufhciently 
small r. 



We now want to get (1721) . Our starting point is Proposition 14.61 and in particular Inequalities 
and (|59[) . The estimates for products and for functions of the density, together with the embedding result 
dSTJ, give us 



W n -\T*) 



T* 



and also 

ll^-l^c^^^ll^llz-,^ 

by use of these two inequalities: 

|| • || L * < || • || fl ._x and HVg^lU- < II^Hb; 

Hence applying (|58|) and (|5^|) to the system ([73]) we deduce 

IKIIz- {B .j < C(||5 n «o|[ fl . ir + Cr) < Cm, 
\\^li UB , r) nL UL ^<Cr + C M C^ 

Hence also ([721 holds true for small r and T*. 
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4.2.2 Step 2 Convergence of the sequence 

In this step we will consider the "difference" sequence 

(5g n ,Su n ,VSn n ) := (g n - g n -\ u n - u n -\ Vtt" - Vtt"- 1 ), Vrc> 1. 
Since taking the difference of the transport term (it + Vfe) • Vw will cause one derivative loss because of 

d 

Vtt" e B^ 1 , we will consider the above difference sequence in the Banach space E* 1 (T*) (recall ([T7]) 
for its definition). 

First of all, by System (J73J), (Sg n , Su n , V6n n ) solves 



(78) 



where 



d t 5g n + u 71 - 1 ■ X78g n - div (« n ~ 1 V5e n ) = F™ -1 , 
a t fa n + (u"" 1 + V6"" 1 ) • V<5u" + \ n ~ l V5ir n = H"- 1 , 
div Su n = 0, 
(S Q n ,Su n )\ t=0 = (A„ eo ,A„ Uo ), 



F"- 1 = -Su 11 - 1 ■ Vg' 1 - 1 + div {{k 71 - 1 - k'^^V^™" 1 ), 

^n-l = h n-l _ _ + y^n-lj . _ ( A «-l _ A „-2) _ 

Next we apply a priori estimates (|48|). ([551) and ([59")) with s = d/p, p G [2,4] and r = 1, to 8g n 
and (Su n , V6w n ) respectively. The use of uniform bounds (fTTj) and (T72"j) for the approximated solutions 
sequence gives us 

(79) H <5fln ll£S?.(fl*^)nLi..(fl*^ +!l ) - C (ll^oH^/P + '^""^^.(s^f)) ' 

(80) ll^llz- (Bp ^ 3 + l|V^l L i t(< /p ni2) < ^([[A^oll^ + \\H n -XuiB^))- 

Now we use the following fact, coming from Proposition 13.101 for any smooth function / = f(p) and a 
given sequence (p m ) m , if we denote f m :— f(p m ) and Sf m := f m — f 771-1 , then 

\\Sf m \\ B ^ < CQ\ e m \\ B . tl ,\\ ( r- 1 \\ B . tl ) \\Sg m \\ B ; tl , if a > ~. 



Therefore, one easily gets 

||*u T, - 1 || fl */»l|Ve n - 1 |L-/, + ||<5e n - 1 |L d /,||e n - 1 |L«./^0, 

Q P.l P.l P.l P.l 

ll^ _1 |liJ,(i#*nLs») ^ C '(II^ 1 |Il^(b^+ 2 ) 

+ f T II^-^I^UVtt"- 1 !!^ + ||<y« n - 1 || fl -/,(||e n - 1 || B -/^ + \\u n -% d/P+1 )). 

Jq p.l P.l P.l P.l P.l 

Hence putting the uniform estimate (fTTj) into (|T9")) entails 

ll*e n || £ oo (s d /»)nLi(B d /»+ 2 ) <C \\\A n g \\ Bd/p +T\\8g n-1 \\ L2 fJ} d/, +1 * + [ WSu 71 - 1 ]] Bd/p \\V g 71 ' 1 ]] Bd/p 
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Apply the above estimate to the term \\6g n II rl , R d/ P +2, appearing in the bound for H n . This makes 
(|5D|) become 

ll^llz^ (B d /p) + ll v,57r "llz^»(B £!/ 1 p ) - C(ll(A n uo, A„_ieo)|| B <i/p +r||(<5e™ _1 ,(5e" _2 )|| L ^^ B d/p + i ) 

+ f \\(6g n ~\5u n -\6u n - 2 )\\ B 4 t , + \\(g n -\g n ' 2 )\\ B ^ P+ 2 + Wu^W^+A. 

Jq PA V Pi 1 Pi 1 P^ 1 / 

Let us now define 

B n (t) := \\Sg n \\ d/p + \\Sg n \\ , d/p+2) + \\Su n \\ d/p + \\V5ir n \\ ! d/p ; 

t V p , 1 / t \ p, 1 ' t V p , 1 / £ v p, 1 / 

then, from previous inequalities we gather 

B n {t) < C IKAn-ifio.Ango.AnUoJILd/, + 

+ r l ( J B"- 1 (t) + S"- 2 (t)) + C / (B 11 - 1 +B n - 2 ) D{a)da, 

Jo 

with ||Z?(t) ||ii([o,T*]) < C- Let us note that, by spectral localization, there exists a constant C > for 
which, for all n > 0, we have 

IKA^cA^o)!^ < C2"Wrt ||(A„eo,A„ Uo )|| L p ■ 

p,i 

Keeping in mind this fact, we claim that the previous estimate implies Y] B n (t) < +oo uniformly in 
[0, T*}. As a matter of fact, for all N > 3 we have 

AT N t N 

Y,B n (t) < C 11(^,^0)11^/, + 2T 1 / 2 Y,B n (t) + / Y,B n (a)D(a)da + B, 

n=l P ' n=l ^° n=l 

where we have denoted by B a constant which depends only on B 1 and B 2 . We can suppose r to have 
be choosen small enough such that, moreover, we can absorb the second term of the right hand side into 
the left hand side. Therefore, Gronwall's inequality entails 

AT 

< G T* ll(eo,«o)|| B -/p , 

n—1 

and passing to the limit for N — > +oo we get our claim. 

Hence, we gather that the sequence (g n ,u n ,Vir n ) is a Cauchy sequence in the functional space 
converges to some (g, u, Vtt), which actually belongs to the space E£ r (T*) by Fatou 
property . Hence, by interpolation, we discover that convergence holds true in every intermediate space 
between r (T*) and E^f{T*) 1 and this is enough to pass to the limit in our equations. So, (g,u, V7r) 
is actually a solution of System ([7]). 



4.3 Uniqueness 

Let us now prove the uniqueness part in Theorem 12. II 

We take two solutions (gi,u±, V7ri) and (ftjj ^2, V^) m Ep r (T*) with the initial data (gi r o,Uifi) for 
i = 1, 2. Then the difference (Sg, Su, \7Sir) = (gi — g2, U\ — u%, V7Ti — V^) solves 



(81) 



dfSg + u% ■ VSg — div (niVSg) = ~Su ■ V g 2 + div ((ki — k 2 )\/ g 2 ), 
d t Su + (ui + V6i) • \7Su + AiV(5tt = hi - h 2 - {Su + V(&i - b 2 )) ■ Vti 2 - (Ai - A 2 )V7r 2 , 

div Su = 0, 

(Sg, Su)\ t= o = (Sg a , Su a ) = (g lfi - g 2 ,o , "1,0 - u 2 ,o) , 
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with the notation = K{gi) and analogous for hi, A; and hi- 

As (gi, ui, V7Ti), (q2, U2, V7T2) G r (T'*), then for any e > to be determined later, there exists T e 
such that 

" ft " i T e (< / i 2+2 )n%(< / 1 p+3 ) ' ll Ui Hi^(B^ +I )' H V7r4 llL^(B^ 1 p+1 ) ^ £ > * = 1,2. 
Let us define (as done before) 

B{t) := ¥Q\\ Lr{B i/ ?) + 11^11^(^+2) + ll*«ll ir(B */P) + \\™*\\ Li (B%?nV>) ' 
Then, the same proof as in paragraph 14.2.21 shows that 

B(t) < C (\\5g \\ B ^ P + \\5uo\\ B 4/,+eB{t)+ J B( s )D(s)ds^j , 

which implies, if e is taken small enough, 

B(t) < e c (\\Sg \\ B d/p + \\5u \\ B d/pj . 

Hence, uniqueness holds for small time. Now, standard continuity arguments show uniqueness on the 
whole time interval [0,T*]. 

5 Proof of Theorems [231 and [27T] 

In this section we try to get a continuation criterion and a lower bound of the lifespan for the local-in-time 
solutions given by Theorem l2.ll It is only a matter of repeating a priori estimates established previously, 
but in an "accurate" way (we use L°°-norm instead of Bp~ -norm) for obtaining the continuation criterion, 
whereas in a "rough" way (we use ([55)) below) for bounding the lifespan from below. 

5.1 Proof of the continuation criterion 

From the arguments after Theorem 12.31 it rests us to prove Lemma 12.41 As argued before Theorem 
12.31 in order to ensure u G L°°(B^ r ), the source terms such as u ■ V 2 b, div (V6 <g) V6) require at least 
u, Vfe G L 00 ^ 00 ) since V 2 6 G ^(B^) only, and also V 2 6 G L 2 (L°°) since Vb G L 2 (B° r ). Similarly, the 
transport term V6 • Vu requires Vu G L 2 (L°°). Therefore in the following, we do not have to optimize 
all the inequalities for the commutators or products and sometimes it simplifies the argument. 



Proof of Lemma \2.4\ It's only a matter of repeating a priori estimates previously established, but in a 
more accurate way. Roughly speaking, in the estimates we use L°°-norm instead of Besov norm, in order 
to require lower regularity. 

Let us consider the density term. Our starting point is (|50p . with this time / = 0, and we argue 
as in proving Proposition 14.51 but control commutators TZj and 1Z 2 (see ((50)) for definition) by use of 
Commutator Estimates (|4"3")) and (|46l) instead. Inequality P"3")) gives us for s > 0, 

|V«||l- || e || B . r + ||Ve||i- ||u|| B . ir ) dr. 

We apply instead (14^1) to control the second commutator term: with a = s, 6 = 1/2 and e > to be fixed 
later, it entails for s > 0, 









/' 


2 js W^W 

L II 'S II LP 




/o 





2 JS 



f \\K 2 \\ LP dr<-[ ||Ve||l„||e|| B . r d7- + e||e||£i (B 
JO f e Jo 



+ 2n 
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So, if e is small enough, p 6 (l,oo), putting these inequalities in (|52p and keeping in mind (1531) . for all 
t £ [0,T[ we get the following estimate: 

( 82 ) IMIZr(Bp,,)nZj(B=+ 2 ) < c(\\g \\ BSpr + \\Q\\ Blr dT + \\V g \\ L ~\\u\\ B s r dT + 



\Vu\\ 



liv e ||5 



Let us now consider velocity field and pressure term: we have 



(83) IMIz r (B- r ) ^ C (\\ U 0\\B lr + \\^\\ZUB }ir ) + \\ h hl(B kr ) + I \\K 3 (T)\\ LP dT 

where we have set, IZj := [u + V6(p) , Aj] ■ Vit. To control the commutator term, a direct application 
of Lemma 14. II yields, for s > 0, 



V s f || [^A^-VuW^dr 
Jo 



< C / \\Vu\\ L ~\\u\\ B s r dT. 

tr JO 



For the second part [Vfe(p) , Aj] ■ Vit, as in the density case, one can resort to Lemma [4.21 to get, for 
s >0, 



V s f \\[Wb(p),Aj]-V U \\ LP dT 
Jo 



< C I ||V 2 p|| ioo \\u\\ B ,dr 



(' I l|Vu||i-||e||fl. ir dr + e || e || E j (B .+ a) 



Let us immediately consider the pressure term: 



l|A(/>)V7r||r ? 



< C\\Vn\\n 



h l(B s r) + ||(A(p)-A(l))V7r|| E j (B 



IVttI 



To bound ||V7r||£i^ ss y we argue as in Subsection 14.1.21 we take the divergence of second equation of 

System ([7)1. we localize in frequencies by operator Aj and we perform a weighed summation. Hence we 
discover, instead of (p5"T) . that 



,u m\ { B^) 



f HA-iVTrllr, + \\h\\ Xl , B . j + ||div((« + V6)-V 

2^ fwdiv^Aj^W^dr 
Jo 



Inequality (|43[) in Lemma |4. 11 entails the control for the commutator term: 



/ HdivaA.AjlVTrJIIipdr 



< 



llVfflli. ||V7r|| B . 



l|Ve|| B .-i HVttH^oo ) rfT". 



By view of (J37J, interpolation inequality helps us to control the above first term on the right hand side 
as follows: 

J HVellio. ||V7r|| B .-idr< ||Ve|U r( iac ) ||V7r|| £ j (B; - 1+ . )) for some fixed ee(0,l), 



<Ce\\Vg\\ 



(s+v)/(l-e). 



£ ll V7r llzi( B|i . 
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with any positive (large) a. We point out here that for any a £ K, there exists C depending on a such 
that the following holds: 

l|A-iV7r|| Ltl(LP) < C||V7r|| £Ks - So) . 

Next, by the divergence-free condition over u and Lemma 14.31 we infer (div u = implies that we 
only need s > 0) 

||div ((« + V6) • V«)|| £tl(B j- rl) < ||V« : V«|| ZKB j- rl) + ||V 2 6 : Vu^^.^ 

< c( I (HVulU- + ||V 2 e|| ioo ) ||«|| B . >p dr + 



o 



t 

2 







The control of the non-linear term h is quite similar as above. We come back to definition (|SJ) and 
consider its terms one by one. As usual, estimation for products in Lemma 14.31 ensures that 



\\ u - y2b hl(B; r ) ~ J ll V2 e|| L oo \\ u \\B i:r dT + \\u\\ L ~(L<*>) 
\(u-VX)Va\\ iUBSp r) < /Vl| 2 L~l|Vg||ic||g|| s , r dr+||g|| £tl(B ,+ 2) + / || Vf?^ ||«|| B . >r dr, 



and 

||(V6-VA)Va|| £ j (B . j < f ||Ve||t-.||e||B. >r d7- + ||e|[£i (B; +> ) 

•J 

Moreover, the last element of h can be treated in the following way: 



||div (V6<8> V6)|| £ i (B . r) = ||A6V6+V6-V 2 6|| £;(Ba } 

f* 

< / ||v 2 e|| LOO \\g\\ B}>r dr + (i + HVelU-^)) llellri^) ■ 

«/ 

Let us collect all these informations: up to multiplication by a constant, we gather 
( 84 ) W^hHB^) £ C(l[V^|U r(£ oo ))S , ( r)||V7r|[ itl(i3 _ So) 

+ / (||Vu||i- + || V 2 g\\ L ^ + \\Vg\\U + HVelltoo) \\u\\ B} dr 
Jo 

+ I (llv w ||i^ + ||u|| 2 00 ||V£ ) || 2 =0 + ||v 2 e || 2 ^ + ||v^|| L ^)||g|| B ^dT 

Jo 

+ (1 + ||Ve|U r(i ») + H|l-(L~)) Ibllzj(B^) > 

and the same control actually holds true also for ||/i||ji( BS y 

In the end, we discover from that ||w|| Zoo ^ BS j satisfies also Inequality (|54"1) . just with an additional 
term ||uo||b« r on the right-hand side. Recalling Estimate (|52l for the density, we replace ||g|| L i^ B s+2^ in 
Inequality (|84[l by the right-hand side of it . 

Thus, we can sum up (1521 and the (modified) estimate for the velocity u, yielding the thesis by 
Gronwall's Lemma. □ 

5.2 Lower bounds for the lifespan of the solution 

The aim of the present subsection is analyzing the lifespan of the solutions to system ([7]). We want to 
show, as carefully as possible, the dependence of the lifespan T on the initial data. This can be done by 
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repeating the a priori estimates previously established, but in a "rough" way. For example, we will use 
thoroughly the following inequalities: 



a2 ||s= . < ||a||i~||a|U5 



(85) \\ab\\ B ^<\\a\\ B; -4b\\ B ^ 
with i = 0, 1, and thanks to Conditions (fT4|) and (Qj) 
(86) 



< 



Mb- 



\\ab\\ L * < \\a\\ L 4b\\ L , < \\a\\ 



Let us point out that, thanks to Theorem 12.31 and embedding results, without any loss of generality 
from now on throughout this subsection we assume 



1 + 



p = 4 



and 



In particular, for any a e R we have || • ||ji^ BCT ) ~ II' \\l\(B" an d this really simplifies our computations. 
We also assume that C always denotes some large enough constant. 

For notational convenience, we define Rq :— ||^o||s= ± ail( l Uq :— ||ixo||b s j, 



From (JH2), we infer that 



and 



U{t) 



Ml 



R(t) + S(t) <c[n a +J^ R(t)(U(t) + 1) dr + R 3 (t) dr^ 



Now, if we define 

(87) T R 

by Gronwall's Lemma, we get 



sup ji > J R 3 (r)dT < 2R \ , 



R(t) + S(t) < CR exp (c(t + J^ U(r)dT 



So, from now on we work with t £ [0, Tr]. 

Let us now focus on the velocity field and the pressure term. We immediately point out that the 
nonlinear term h will make density terms with critical regularity appear. Therefore, unlike what has 
been done in previous sections, we decided not to use systematically interpolation inequalities to isolate 
the term S(t): we consider instead S'(t), which controls high regularity of the density. 

Our starting point is estimate 



AVttI 



L *( B li) 



2-'' 



f\\Ki{r)\\Lv 
Jo 



dr 



where, as before, IZj := [u + Vb(p), Aj] ■ Vu. Now we apply the classical results, e.g. lemma 2.100 of 
[?] , to control the commutator term, and we get 

(89) U(t) < c(u Q + \\h\\ B ; A dT + ||AV7r|| B?i dr + J U 2 (r)dT + V "(t) S' '(t) dr^j . 

Let us now establish some stationary estimates, which will be useful to complete the bound for U(t). 
We will often use (|85l) to deal with products. 
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(i) Let us start with r - We consider its terms one by one, and we get 

||u-V 2 6|L < CUS' 
|| (it- VA)Va|| BS < CURS' 



p,i 



||(V6-VA)Va|| B . < CR 2 S' 
||div (V6® V6)|| B . i < 2 ||V 2 6® V6|| B3 <CRS'. 
Therefore, one infers that 

\\h\\ B;i < C(US' + URS' + R 2 S' + RS') . 

(ii) Now we focus on the Besov norm of the pressure term. First of all, we have 

||AV7r|| B?ii < C(l + iJ)||V7r|| fl . ir . 

Now, we can handle ||V7t||b« r . Since r = 1 here, we do not have to worry about the order of taking 
the integration with respect to time L\ and taking the l r norm of the dyadic sequence. 

Hence, we can first "forget about" the time variable and discuss ||V7r||.g* . In order to control it, 
we separate low and high frequencies and use Bernstein's inequalities: we get 



IIVttIIb^ < HA-xVttH^ + || (Id - A^)\7n\\ B s p i < C [\\Vn\\ L 2 + UAttU^-x 

The equation for the Laplacian of the pressure comes from (|64[) : 

-Att = V(logA) ■ Vtt + A _1 div (v ■ Vit - h) , 
from which it immediately follows 

IIATrll^-a < ||V(logA) • Vtt|| b . +(y| fl --i+l)(||div(i;-V«)|| fl .-i + ||/i|| B;i ) 

< \\Vq\\ b s-s ||V7r|| B .-i + (||e|| fl .-i + 1) (||Vt> : V«|| B .-i + \\h\\ B . ) 

p,i p,i p,i \ p,i p>^/ 

(where we used also that divit = 0). So, applying the interpolation inequality for the embeddings 
r B^ 1 B p then the inclusion L 2 B p ^ 2 and finally the Young inequality, we gather 

||V7r|| B . t < C ((1 + R s ) ||V7r|| L2 + (1 + R)\\Vv : V«|| fl .-i + (1 + R)\\h\\ B;i ) , 
for some 8 > 1, thanks to divit = 0. Then, we have also 

(90) ||AV7r|| B|il < c[(l + R)(l + R s ) ||V7r|| ia + (1 + R) 2 U(U + S') + (1 + R) 2 \\h\\ B; J . 

(iii) Now we focus on the L 2 norm of the pressure term. Applying Lemma 13.131 to equation (I64p 
immediately gives 

IIVttH^ < C(\\h\\z.» + ||fVu|| L2 ) . 

Thanks to ([55]). we find 

||u-Vu||x,2 < C\\u + V6|| BS -i ||Vit|| B a-i <C(U 2 + UR) . 

Using also (|85[l and (|55|) . it's easy to control the L 2 norm of h: 

\\u-V 2 b\\ L2 < C \\u\\ BS -i \\V 2 b\\ B s-i < CU S' 
||( M -VA)Va|| L2 < CH| B .-x||Ve||^_x < CUR 2 
||(V6-VA)Va|| i2 < C\\Vg\\l^ < C R 3 

PA 

||div (V6® V6)|| ia < 2 || V 2 &<g>V&|L 2 < CR S' . 
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Before putting all these inequalities together into (|89|) . let us note the following point. Due to the 
fact that 6 > 1, for any m > we have 

R m (1 + R s ) <C m (l+ R m+S ) , 

and we have to deal with only a finite number of powers of R (the biggest m is actually 4, by previous 
stationary estimates) . Hence, if we define (denoting by r\ some big enough exponent) 



£{t) := 1 + Rl cxp^C* (t + J o U(t)(1t 



thanks to (155)) we can bound all the terms of the form R m and R m (l + R s ) which occurs in our estimate 
by C £ , with C > suitably large. 

Now, we are ready to complete the bound for ||A V7t||bs i . From the previous step (hi), we get 

(1 + R){1 + R 5 ) \\Vw\\ L 2 <C(£U 2 +£U + £US' + £ S') , 
while from the step (i) we infer 

{l + R 2 ) \\h\\ B ' ir < C(£US' + £ S'). 

Then from (|90l) we finally get 

IIAVvtIIbs < CI £U 2 + £U + £U S' + £ S' 



Let us note that, as £ > 1, we can bound also by the previous quantity. 

Now we put this last inequality in (|89p . Using again the fact that £ > 1, we find 

(91) U(t) < C (u + J (£U 2 + £U + £US' + £ S') dr 
for all t G [0,Tr]. Therefore, if we define 

(92) T v := sup 1 1 > 



(£ U 2 + £ U + £ U S' + £ S') dr < 2U 

then in [0,Tr] (~1 [0,Tjj] we get U(t) < CUq and we manage to close the estimates. 

So, our next goal is to prove that, if we define T as the quantity in (|23p . then T < min {Tr, Tjj}, i.e. 
both conditions in ([87)1 and ([9"2")) are fulfilled. 

Let's first tackle the case Uo = 1 and then we will see how to deal with the general case by use of 
Proposition H31 First of all, from ([55)1. ([9"T)l and the definition of £, in the interval [0,T R ] n [0,Tu] we 
have (if Uq = I) 

(93) R(t) + S(t) < CR e ct and £{t) < (l + R n Q ) e ct , 

for some suitable positive constant C. Note that in [0, Tr] D [0, Tjj] we have 

/ (£U 2 + £ U + £ US' + £ S')dr < 2C (1 + R% +1 ) e 2Ct . 
Jo 

Therefore, Tu, defined by ([9^| . is bigger than any time t for which the quantity 2C (1 + Rq +1 ) e 2Ct is 
controlled by 2. Hence, for suitable values of L and with 3t = r\ + 1, Uq = 1 in ([2~3")l . we have that T < Tjj. 
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Now, let us consider Tr. By (|93| . it's easy to see that the left-hand side of the condition in (|87| can be 
controlled in the following way: 



f R 3 ( T )dT < CRl f e 3CW dT < c|| 



o e 3Ct 



Hence, Tr is greater than any time t for which Rq e 3Mt < C , for some convenient constant C. But this 
condition is always verified in [0, T], up to change the values of L in (|2"5|) (recall that X > 2). Therefore, 
in the end we gather that T < min {Tr, Tjj}, and Theorem 12.71 is proved for Uq = 1. 

For arbitrary initial data po ~ b u o hi -E^ - ^ 4 , we have a unique local-in-time solution (p, u, \7ir). If 
we now set 

= ll u o|| B i+<i/4 , 

4,1 

then we can see (p, u, V7r) as a rescaling of some solution (p, u, Vtt) having initial data 

(pq,u ^(x) := ^p ,£ _1 wo) (e _1 a;) ■ 
In particular C/o := ||uo|| R i+d/4 = 1 and hence the lifespan of the solution Co, u, V7r) is bigger than 

4,1 



ilog (TTlP TP ) =il °g(- 

V1+ \\po v x 



/4 y v l + (e 1 ||po - l\\ B i+d/i)~ 

In virtue of Proposition 12.51 the lifespan of the solution (p, u, Vir) is larger than 

L f L 
T^log 



n i+(^ 1 ||P0-l]] Bl+ ./4) 3£ 



This completes the proof of Theorem 12.71 thanks to Remark 12.81 



6 Other cases 

We want now to deal with finite energy initial data and with small initial densities. Before passing to 
the proof of Theorems 12.91 and 12.121 we need an estimate for parabolic equations in Besov spaces of type 
r , in order to tackle also the case p — oo. 

6.1 Linearized parabolic equations in B s oor 

The aim of this paragraph is to show an estimate for parabolic equations in Besov spaces of type r . 
Such an estimate will be fundamental in the sequel. 

In the present subsection we resort to notations and tools from Homogeneous Paradifferential Calculus, 
which enjoy most of the properties we have seen for the non- homogeneous case (see also Section [3]) . Wc 
refer to Chapter 2 of [?] for a detailed description. 

The result we want to prove in the present subsection is the following. 



Proposition 6.1. Let p 6 S solve the following linear parabolic-type system 
(94) 

with k, f, po G S and 



d t p - div (kV/j) = /, 
p\t=o = Po, 



< p* < po < p* , < < n(t, x) < k* 
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If s > 0, r 6 [1, +oo], then the following estimate holds true: 

(95) IIPlll r(fl|o , r )nZj( B S+?) < Cz CK ~ {t) >< (lIPolUso,, + H/U^J . 

where C is a constant depending on d, s,r, p*, p* , k* and ||rt||^c»(c;«) / or some e £ (0,1), and 

#«,(*) := /' (l + [|V«|[i„ + || V K ||™ x f /(1+s) ' 1} ) dr . 

Remark 6.2. Let us point out that the above proposition provides us with estimates analogous to those 
in Proposition 14.51 (for instance, one can compare the expression (|4T)|) with the definition of -Kqo)- 
Although here we are considering a parabolic equation with no transport term, the result can be eas- 
ily extended to that case by a direct use of product estimates as well as interpolation and Gronwall 
inequalities. 

We will first focus on the case s S (0, 1) and r = oo, i.e. we take solutions in time-dependent Holder's 
spaces L^(C e ) Pi L^(C 2+C ), e e (0, 1). Note that Maximum Principle applied to parabolic equations 
has already given us the control on low frequencies of the solution: 




besides, the classical a priori estimates for heat equations are simpler (at least formally, see Proposition 
16. 3p in homogeneous Besov spaces. For these reasons, we only have to focus on the homogeneous Holder's 
space: 

E e := Zf?(C e ) nZ^((7 2+£ ), £ e(0,l). 

We can localize the function p into countable functions g n , each of which is supported on some ball 
B(x ni 5), with the small radius <5 6 (0, 1) to be dertermined in the proof. It is useful to notice that time- 
dependent Holder's spaces may be described in terms of finite differences (see Proposition 16.51 below). 
Roughly speaking, up to a constant, the E^-norm can be determined locally, which reduces to bound 
{g n } instead of the whole p. The systems for {g n } (see System (|101l) below) derive from (|94p after 
mutiplying by some partition of unity and thus the coefficient k can be viewed as a small perturbation 
of a function which depends only on time t, because g n is supported on a ball with sufficiently small 
radius. Consequently, changing the time variable and making use of estimates for heat equations entail 
an a priori estimate for g n in E e , with the "source" terms being either small or of lower regularity, and 
hence easy to control. At last we will show how the Holder case will yield the general one. 

We agree that in this subsection {g n (t,x)} always denote localized functions of p(t,x) in x-space, 
while pj as usual, denotes Ajp (localization in the phase space). 

Before going on, we recall a priori estimates in homogeneous Holder's spaces for heat equations with 
constant coefficients (see Chapter 3 of [?]): 

Proposition 6.3. For any s 6 K, there exists a constant Co such that 
(97) \\n L ~ ( c*)nZUC*+>) ^ Codl/ollc, + ||/|| Z x (<>) ), 

where fo,f,F e t S(M <1 ) and are linked by the relation 

F(t,x)=e tA f Q + [ e^ A f(r)dr. 
Jo 

It is also convenient to show an a priori estimate for the paraproduct T v u in the space L)p(C s ), similar 
to Lemma POl 
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Lemma 6.4. For any s > 0, e > 0, a, b > 0, there exists a constant C £ ~ e a ^ b such that 



(98) \\ T M\l U c>) < c cJ q IK*)II£ \Mc-°+e\Mz U c>>y 

Proof. First, let us notice that, 



\\T v u\\ Zu6e) < 



[ ^IHIc»ll A ^lk° 



For any e > 0, a, b > and t <E (0, T), we fix an integer 



N t = 



iiog^Ki)!!^) 



then noticing that e 1 ||u(t)|| ( j s ~ 2 * , we have 

i-T r T 



f ^WnW^WAjvU- < j Y, 2 ia |Kt)||c s 2-^||A i «|| i o + 53 2-^||«(t)|| 6 .^ 6 ||A i «|| £ o 
< f 2 N * a \\u{t)\\cM\c-« + l- {j - Nt)h eV h \\^v\\L°° 

J ° j>N t + l 

Jo ^ 

f T \\u(t)\\^ /b \\v\\ d . a +esu V fy b \\k 3 v\\ L ~. 
Jo j Jo 



Thus the lemma follows. 



□ 



Recall now the characterization of homogeneous Holder spaces and time-dependent homogeneous 
Holder spaces: 

Proposition 6.5. Ve G (0, 1), there exists a constant C such that for all u G S, 

\\u(x + y) -u(x)\\ L o 



(99) 

and 
(100) 



C- X \\u\\c< < 



c-'IMIzjcd.) < 



\u(t, x + y) — u(t, x) 



<C\\u\\ dc , 



dr 



<C\\u\\ Z} 



Proof. The proof of can be found in Chap. 2 of [?]. 
Let us just show the left-hand inequality of (11001) . Since 



A jU (t,x) = 2 jd / h{2 j y)(u(t,x-y)-u(t,x))dy, 



then 



V[|A i «(T 1 .)|U r <2*' / 2*\ V \<\hplv)\ [* MT ' X y [ Ae U{T ' X)h * drdy 



" 



\y\< 



< c 



1 \\u(t, x + y) — u(t, x) 



which ensures the left -hand inequality of (|100[) . The inverse inequality follows immediately after similar 
changes with respect to time in the classical proof. □ 
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Next we will prove Proposition 16.11 in three steps. In Step 1 we will mainly deal with the localized 
solutions {g n } or equivalently, after a one-to-one transformation in time variables, with {g n }, which solve 
heat equations. Thus Proposition 16.31 ensures an a priori estimate for g„ and Lemma 16.41 provides the 
control on "source" terms. Finally, thanks to Proposition 16.51 we can carry the results from {g n } to p: 
this will be done in Step 2. In order to handle general Besov spaces of form B^ X}r , we again localize 
the system, but in Fourier variables, in Step 3. Then we apply the result of Step 2 to pj and a careful 
calculation on commutator terms will yield the thesis. 

Step 1 Estimate for g n in E e 

Let us take first a smooth (e.g. C 3 is enough) partition of unity {^ n }neN subordinated to a locally finite 
covering of R d . We suppose also that the i[i n 's satisfy the following conditions: 

(i) Supp^n C B(x n ,5) = B ni Vn G N, with i5 < 1 to be determined later; 

(ii) En^n = 1; 

(iii) < ip n < 1, Vn S N with ij} n = 1 on B(x n ,S/2); 

(iv) ||V^n|U~ < C|*|-I"l, Vn G N, for |r?| < 3; 

(v) for each x G M. d , there are Nd (depending on the dimension d) elements in {ipn}nefi covering the 
ball B(x,S/2). 

Now by multiplying ip n to Equation (|94p . we get the equation for compactly supported function 
g n = pip n which is supported on B n : 



(101) 
where 



dtQn - n n Ag n = (k - K n )Ag n + Vn ■ V 'g„ 

g n \t=0 = Q0,n = IpnPO, 



vol(±l n ) J Bn 



is a function depending only on t, and 

(102) g n = -2kVV>„ • Vp - (nAipn + Vk • Vip n )p + fip n . 

For convenience we suppose that there exists a positive constant C K ~ I k IIl°°(<7c) such that 

(103) \n{t,x) - K (t,y)\ < C K \x-y\\ Var, y G M d , t G [0, to]. 
Noticing that, by (|103[) . we have R n > 0, this ensures that, for all t G [0,io], 

1 



(104) ||«/k„-1|U»(b„) <k* 



ABn) 



K(t,x) - K{t,y)dy 



L°°(B„) 



In order to get rid of the variable coefficient K n (t), let us make the one-to-one change in time variable 
(105) T^r{t) = f K n {t')dt'. 

Therefore, the new unknown 



g n (T,x) = p n (t,x), 



satisfies (observe that ^ = n n {t)) 

(106) { °re» A ^ = {Sh- 1 ) A ?« + W\ ■ + Wr 

y gn\r=0 = P0,n, 

where k(t,x) = K(t,x), k„(t) = fi n (t), p(r,x) = p(t,x), g n (T,x) = g„(t,x). 
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This is a heat equation and hence by view of Proposition 16.31 it rests to bound the "source" terms. 
Estimate (j9"5)) and the following one, 

\\f u V + R(u,v)\\ i}r{Gc) < C||u|| £ o ?(2 ;oo ) ||«|| £ i (<: ,. )1 

imply that the first source term of Equation (|106[) can be controlled by 

"k(t,-) 



k„(t) 



- 1 Ag n (r, 



< C \\K/K n ~ l|lioo(i«.(B„)) ll^5r»Hii 

+ C m / H/s/rtn - ||Ap„|| d , e _ 2 +r?i||A£> n || £ i i(6e) , 



for any r/i G (0,1) with C m ~ ?y x E . Besides, Inequality (|104[) ensures that for all t G [0,to], with 
to = r(t ), 

\\n/K n - l[[i«(B n ) < C k k7M% 

which implies, for T £ [0,to], 
'k(t,-) 



(107) 



k„(t) 



- 1 Ae„(r, 



<C m || «|U e Iknl! c< + (CC K K, <5 £ +?7l)||£> n || £ i ( p 2+£) . 



By Lemma POl for any 77 > 0, there exists C, ~ rj 1 such that 



||T„u + i?(M,w)|| z ^ (d , e) < C*^ y ||w||i«>||t;||c«-i + 7 ?b|| £ i, ( c;e+i)- 

Thus, also by use of Lemma 16.41 with a = 1 — e and 6 = 1 + e, for any 772 £ (0, 1) we have the following 
(with a, 2 - r^ 1 ): 



(108) 



Vk(t, •) 



V£>n(r, •) 



«n(r) 

Now let us choose 5, r\\ , 772 such that 



< a )2 / HV/ellioo + ||Vk||^ + ijallft.lb (62+e) . 

ij.(6«) «/o v ; ' 



(109) 



CqCC^S 1 , Cqt?! , C 77 2 < 1/6, 



with the same Co in (|9"T|) . Then, from Proposition 16.31 and estimates (|107[) . (|108p . for any t G [0,Xb] we 
get, for some "harmless" constant still denoted by C, 



IKIIz~(Ce)nZ^+e) ^ C ^lleCnllo- + ^ (JMI^ + ll VK lli~ + H Vk II£* J llf?«llc« + llSnllz^) 

Since k, < k„ < k*, after transformation in time (|105[) we arrive at 



(110) 



\\e™\\L™(C')nL 1 T (C 2 +') - C \\80,n\\d* + / ■ K l||0n|lc« + 



for all T G [0,t ], with 



' x > = ll«ll^ + llv«nio + ||v«||vi* ! 



provided 

(111) 6~ e — C'sC K < C5 1| K || L oo for some constant C$ depending only on d, e. 
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Step 2 Estimate for p 

Now we come back to consider p = g n . By assumptions on the partition of unity {VVi}, for any x 
there exist Nd balls of our covering which cover the small ball B(x, 5/ A). Therefore, from Inequality (|M|) 
we have 



\\p(r,x + y) - p(r,x) 



\y\< 



< C sup / — + C sup / j-rj - 

\v\>s/iJo W \y\<s/iJo \v\ 

whose second term can be controlled by 

rt su Pn \\Qn{x + y)- Qn{x + z)\\ L «= 



N d C sup 

\y-z\<S/AJQ 



\y-zY 



Thus we find 

IMIzi ( ce) - C6 C J Q WpW^ + N dCsnp\\g n \\ llt{dt) . 

Similarly, we have 

||HIZ|o (c?e) < C , <T e ||p|U t ~(z,«>) +Csup||e n |ll«» ( c.e). 
Since V 2 p = X)n(^ 2 ft 1 )> f rom the same arguments as before we infer 

HHb(c^) ^ cII^pIIzhc') ^ / llv 2 HU~ + c sup ll^ki,^. 

Therefore, to sum up, for all t 6 [0,io], 



<cr f (||p||L r ( L ~) + ^||v 2 plU») + 
+ Csup Mleo.nllce +^ Kx\\Q n \\ de + \\g n \\zi(C')j > 

with the second inequality deriving from Estimate (jllOl) . Thanks to ((96)1 and the fact that 

llft.Ho- = |Hn||c < C[[p|| C «|[Vn||c« < C5^\\p\\ c ,, 

we thus have the following estimate for p in the nonhomogeneous Holder space 



ll(c 2 +>) - c \\p\\l^(l°°) + C J 



< cr £ M|p ||c* + y o (liv^iu- + \\ P \\ L ~ + n/|| L oo) + ^ ^hhIc- 

+ Csup 

It rests us to bound g n uniformly. In fact, starting from definition (|102[) of g n , we follow the same 
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method to get (|107l) and (I108P to arrive at 

Il0nllri ( c.) < c J o ^"Hllwv^iil- + Hkv^h^) + ||kA^„ + vk • wgu«^ \\ P \\ d . 

+ v\\p\\ 

< cv~ 2 jf (l + ||«||^ + ||«|| C i+.) IMIc + v\\p\\zi { c^) + C6-<\\f\\ Znc .y 

where we have used \\f\\ L i (La o )n Zi^) ~ C \\fh\(f!<>y 
We finally get a priori estimates for p: 

IIPllz t -(c«)nZj(c^«) ^ CS ~' (\\P°Wc- + J \\P\\C + WfWlKcn) + C5 ~ 2 j K *\\P\\C*, 

with 

K 2 = l + \\ K \\& m >c(k 1 + 1+ \\k\\^ + ||«|| c i+. 
Thus, by a direct interpolation inequality, that is to say 

8~ e \\p\\Ll(C*) < C v S ~ 2 ( WpWc' + v\\p\\l\{c^+'V 



(i 



Gronwall's Inequality tells us 

(112) Wphnc^ZKc^) < CT " £ex p ( cs ~ 3 f a K *) + Wfhuo) > 

which is just Estimate (|95|1 when s = e and r = oo. 
Step 3 The general case B^ r 

Now we want to deal with the general case B^ r . Let us apply Aj — Aj-i + Aj + Aj + i , j > 0, to System 
(EH), yielding 



(113) 
with 



dt^ - div (kVpj) = fj - Rj, 

Pj\t=0 = Po,j: 



Pj = Ajp, fj = Ajf, Rj = div ([«, Aj]Vp), p 0j = Ajpo- 

We apply the a priori estimate (|112p to the solution p 3 of System (|113[) . for some positive e < min{s, 1}, 
entailing 

WPikriOnZHC*.) ^ CS ~ e ™P [C5- 2 J^K^j (||p 0J || Ce + l^-^-llzj^) . 
Let us notice that for j > 0, denoted by pj = Ajp and p q = A q p as usual, then we have 

Ajpj = pj and A q ~pj = if |g — j\ > 2. 
Hence the above inequality gives 

2 je ||p J IU r( ^) + 2^ 2+£ ) /" Mi- < 



o 



r-2 



(** E (i 

\ b'-?|<i 



<C<r e exp(C<r^ K 2 j 2^ ^ (Hpo.flHi- + y o ||/ 9 |k«)+ 
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Let us focus on H-Rjllji^e) f° r a while. As usual, we decompose Rj into four parts: 
R)= div ([*Vi«>A/]VAg/>), R 2 j= divA5(A 9 KS* g+1 Vp), 

|j-«l<3 1>j-3 

Rj = div (A q K AjSq+i Vp), Rj = div ([A-\K, Aj] Vp), with k = k — A_ik, . 

9>J-3 

1 2 

It is easy to see that the Fourier transform of the terms Rj and JL is supported near the ring 2 3 C; thus, 
by Lemma 2.100 of [?], we get (for some sequence (cj)^ 6 £ r ) 

\\R]\\ Zuce) <CVW [\\[S q -iK,Ai]V\p\\L- 

\j-l\<3 J ° 

<C2>* I ||V«|U-||VA,p|| L - 

\j- q \<z Ja 



<C2^-^Cj [ ||V«||i-||Vp[| %p 
Jo 



and, for some s > — 1 and (cj) . S i? r 



||^ 2 || ZKCe) <C2^ 1+e ) £ Aa^iu-iis^vpIU- 

9>i-3 J ° 

1>j-3 J ° 

[ ||V«|| fli , r ||Vp|| i -. 
Jo 



Let us now consider Rj, Each A q K AjS q +iVp has Fourier transform supported near a ball centered at 
origin with radius of size 2 9 . Therefore, arguing as above, we find, for some (cj)j G £ r 

ll^llri(c)^ E ., su p 2J ' (1+e) f \\& q K\\L~\\AjVp\\ L ~ 

<C 2?(£+1) r2- 9 ||V«|| L oo Cj 2^ s ||Vp|| BSoi . 

<C2^-^Cj [ ||VK||i-.||Vp|| BSer . 
Jo 

Since Fourier transform is always supported near a ring 2-?C, arguing for instance as in the proof 
of Lemma B~2l (see the Appendix), we can bound \\Rj || a ^ so by the arj ove quantity. 

To conclude, we have got a priori estimate for pj\ 

K 2 | • 



||PilU- ( L-)+ 2 2 ^ llp^li- < C<T e exp (CS- 2 J\ 

E / ||/ fl IU-)+2- i 'c j / ||V«|| B!o J|Vp||lo.+ / ||V/s||z-||Vp||j 

. „l^n JO 7 JO JO 



b'-«|<a 
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Therefore, we multiply both sides by 2 JS (for s > —1) and then take l r norm, to arrive at 
IMIz r(S ^nLi(^) <C5-*eiq>(c6- 2 J*K^j x 



x (\\P°\\B L , r + \\f\\z l{BL r) + J* \\VK\\i 



\\Vp\\ B . +||V«|| B . J|Vp|| 



The direct application of interpolation inequalities and embedding results reduces the above estimate 
into pjj) . 

6.2 The case of finite energy data 

The proof of Theorem l2.9l is just as the proof presented in Subsection l4.2l with some changes pertaining to 
the energy. We will follow the standard process in proving the local existence result also: we construct a 
sequence of approximate smooth solutions which have uniform bounds and then we show the convergence 
to a unique solution. In each step we will try to sketch the analysis first and then present the proof into 
details. 

Let us make some simplification. Sometimes a few estimates may depend on the existing time T* , and 
hence a priori we suppose that T* < 1. We also assume that all the constants appearing in the sequel, 
such as C, Cm, Ce, are bigger than 1. We always denote Sb n = b(p n ) — b(p n ~ 1 ) and 5a n = a(p n ) — a( ( o n_1 ). 

6.2.1 Construction of a sequence of approximate solutions 

As usual, after fixing (g° , u°, V7r°) = (go, uq, 0), we consider inductively the n-th approximate density g n 
to be the unique global solution of the following linear system 

, , ( d t g n + u n - 1 ■ \7g n -div(n n - 1 Vg n ) = 0, 

L g |t=o — Qo, 

with k" -1 = n(p n ^ 1 ) and the n-th approximate velocity and pressure (u n , V7r Tl ) satisfying 

( d t u n + (u 71 - 1 - n n - 1 {p n )- 1 Vp n ) ■ Vu" + A"V7r" = h n -\ 

(115) I divu™ = 0, 
[ u™| i= o = wo, 

where we have set as before, A" = X(p n ) and the same h n ~ 1 defined in (1731 : 

(116) h 11 - 1 = ( /9 ,l - 1 )- 1 (A6( j o"- 1 )Va( j o"- 1 ) + li"" 1 • V 2 a(/°" _1 ) + Vfe^"" 1 ) • V 2 a(/>™ _1 )) ■ 

We pay attention that, compared with System ([731 . the coefficients it" -1 — K n ~ 1 (p n )~ 1 Vp n and A™ of 
System (|115[) here are chosen to keep accordance with Equation (|114[) , in order to get the energy identity 
for u n . Indeed, noticing divu™ = and Equation (|114[) for p", we can take the i 2 (R d )-inner product 
between Equation (|115p i and p n u n , getting at least formally □ 

(117) ~ { p"K| 2 = [ p n h n ~ l -u n . 

i at j R d j R d 

Furthermore, the initial data for {g n } are chosen the same. We will see later that this choice ensures us to 
estimate the difference sequence {Sg n — g n — g n ~ 1 } n >2 hi Space Ct* [H 1 ) (T* denotes the existence time), 



2 We can also follow the scheme 1731 which, gives the following by taking L 2 -inner product between 1731 12 and p n 1 u n 

-— [ p »-i|«"| a + i / divfp"- 1 ^™- 2 -u™- 1 ) - (ft"- 2 - K^-^Vp"- 1 )!^! 2 = / p 71 - 1 ^- 1 ■ u n . 

2 dt J R d 2 J R d V J J Rd 

For simplicity we choose iterative linear systems (11141) and I H15I here. 
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although the initial density go belongs to L 2 only. This estimate makes it possible to bound the difference 
of the source term h n — h n ~ x (notice it involves terms such as A5b n \7a") in Space L\* (L 2 ). Therefore the 
convergence of the velocity sequence in Space Ct*(L 2 ) follows. We will explain the convergence process 
in detail in Subsection 16.2.21 below. 

In this paragraph, we aim at proving the existence of the solution sequence (p™, u n , V7r") and uniform 
estimates for it. We want to show estimates (|7T|) and (|72")l . with a change pertaining to Vir n : 

(118) ||V7r"|| £ ^ (B , r)nL ^ (L 2 ) < r^, with s p = maxjs , s - ~ + ~ J . 

We also prove the following inductive estimate involving energy: 

(H9) ||e*|U- (£.) + l|Vg"|| L| , (i2) + ||« n |U~ < c e e 0i 

with Eq := 1 1 go ||i 2 + || u o||l 2 and some constant Ce depending on d, s, p*,p*. 

In fact, the subtlety still comes out when dealing with the pressure term V7T ra . Getting inductive 
estimate (|118|) relies on the divergence- free condition of u n , which helps us to write the equations for 
7T™ in different forms, such as (|121[) . (|122[) and (|123[) in the following. Informations on low and high 
frequencies issue from these equations separately. This yields estimates for V-7T itself finally, thanks to 

(120) HV7r|| Z j (B . r) < C||A_iV7r|| L i (w) + \\ATr\\ ZliB ,- ly 

In Subsection 14.21 where p £ [2,4], by view of the divergence-form elliptic equation for 7r n , the 
quantity ||V7r" || (z, 3 ) (involving the low frequency information) can be controlled by a simple use of 
Holder's Inequality ||/.g||L 2 < ||/||.L4||<7||z,4 to the quadratic "source" terms. Now, we also apply "div" to 
Equation (j 11 5[) . getting the following equation for 7r™: 

(121) div (A"V7r") = div (h n ~ l - (u"- 1 - fs n - 1 (p T, )- 1 Vp n ) • V«j . 

But here, thanks to Estimate (j 1 19[) . one uses ||/g||L 2 < ||/||l 2 II^IU 00 to bound || Vtt™ || x, 1 » (i 2 ) ail( l hence 
the low frequencies, when p > 2. 

If p £ (1,2), instead, we multiply Equation (|115p i by p n and then apply "div" to it. Recalling that 
div9tu n = 0, we get a Laplace equation for 7r Tl : 

(122) Att" = -div ((p n - l)d t u n + {p n v n - 1 - n n - l Vp n ) • Vit") + div (p n h n - 1 ). 
One observes that if p £ (1, 2), then 

B;- l - d / q ^ L q , Vq > 2. 

Thus we perform ||/<7||lp < |/||l 2 H^Ulp* with 1/p* + 1/2 = 1/p to these right-hand side quadratic terms. 
Since g n is already in L q for all q > 2 according to Maximum Principle and Energy Estimate, this requires 
||9tu Tl || i i( i 2). Luckily it is related to || V7r ra || L i( £ 2) by Equation f| 1 15[) x . 

In order to control the high frequency, it is enough to show A-7r n in L\ (B*" 1 ). Let us rewrite Equation 
([T2T|l as 

(123) Att" = Vlogp" • Vtt" + p n div (V 1 " 1 - {u^ 1 - n n ~ x {p n )-^ p n ) • V«j. 

The quantity || V7r n || can be interpolated between L\(B^ r ) and L\{L P ) when p £ (1,2) (resp. 

L\(L 2 ) when p £ [2, oo) since L 2 d ^ 2 ) and the index s p will come out from these interpolation. 

As in [?], let us also notice that the limit functional space with s — r = 1 is no longer an 

algebra as spaces B^. 1 with p < +oo. But we still have the following product estimates by considering 
paraproducts and remainder separately: 

(124) \\f9\\Bl A <C\\f\\B L Jg\\ B ^ Ve>0. 
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Thus we just have to modify the above interpolation a little to make quantity ||V7r"|| ^ appear. 
We refer also to [?] for the inequality (due to div u n — 0) 

(125) l|div(u"- 1 -Vtt n )|| B o oil <C||u B - 1 || B i oil ||Vu B || B o 0ii . 

Now let us realize the above analysis. 
Since (g°, u°, Vtt°) = (g , u , 0), then by choosing small T* , Estimates (J7T]), d]), (fTT5|) and ^TW\i all hold 
for n = 0. Next we suppose (g 11 ^ 1 , u™ -1 , V7r™ -1 ) to belong to the functional space E defined by 

(126) (c(M+;^ J ,nL 2 )nLL( J ff 1 )nZ 1 1 oc (Bp+ 2 )) x (c(R + ;B;, r nL 2 )) x (Z\ oc {B s p , r ) r\L\ oc {L 2 )), 

and such that the inductive assumptions hold. We just have to show that the n-th unknown (g n ,u n , V7r n ) 
defined by System (|114[) and (|115D belongs to the same space, satisfying the same conditions. 

Instead of working with g n and ti" together, we first consider g n independently. By smoothing out 
the initial datum and coefficients, applying Proposition 14.51 or Proposition 16.11 and then showing the 
convergence, we can get the unique global solution g n of the linear system (|114l) . The process is quite 
standard and we omit it. One observes that Estimates £25]) and (JM]) imply g n £ L^°{B s p r ) n Zj(i?j;+ 2 ) 
for any finite t > 0. Thus in particular g n £ C(M. + ;Bp r ). 

On the other hand, since u™~ 1 £ L^ C (L°°) and div it" -1 = (noticing Equation (|115p 2 ). energy 
inequality $Q for g n follows. Thus g n 6 C(R+;L 2 ) n Lf oc (H 1 ), and we have 

\\Q n \\L^) + \\^Q n \\ L l,(L^<\c E \\g Q \\ L ,. 
It is easy to see that, by Maximum Principle, 

P* < 1 + g n {t,x) = P n (t, x) < p*, vt el + ,ie M. d . 

As in Subsection 14.21 we introduce g^ to be the solution of the free heat equation with initial datum 
go, which satisfies (|74l) and (|75|) also. Correspondingly, the remainder g n := g n — g^ solves System (|77p. 
with gL instead of g\. Propositions 14.51 and 16.11 thus imply that for p £ (1, oo], 

where 

/" = -u 11 - 1 ■ Vg n - u n - x ■ Vg L + div ((/t"" 1 - l)Vg L ). 
Here C" _1 (t) depends also on l^™ -1 ) when p — oo, and by embeddings, we can take 

JC n -\t) : =t+\\V K n - 1 \\ 2 LUL~)+ IIv^-^Ukb-j. 
Since the inductive assumption (|71l) holds, we derive on the whole time interval [0, T*] 

C n - 1 (t)e c "" 1 «' cn ~ 1 W <C K , 

for some constant C/c depending onljj^l on M. Furthermore, product estimates, interpolation inequality 
and Estimate (|75|) ensure that 

Therefore the smallness statement (I76p pertaining to g n is verified and hence inductive assumption (I7ip 
holds for g". 

3 In fact, Cic also depends on T* , which can be "omitted" since we have supposed a priori T* < 1. 
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As to solve System (j 1 1 5[1 . a convenient way is to view it as a transport equation of the velocity 
u n . For each finite time t, if there exists a constant Ct (depending on t, ||(£>, it™ , Vir n ~ 1 )\\ Es ( t j, 

W^WL^iL^nL 2 ^ 1 )-! k = n — l,n, and ||u n_1 Hl^l 2 )) such that the L\ (_B| r )-norm of the "source" term 
-A™vV n + is bounded by C t (l + U n {t)) with 



Wit) = \\u n \\ 



then a standard proof gives a unique global solution u n £ C(E + ; r ) (see the proof of Theorem 3.19 
in [?]). If furthermore Vtt™, /i™ £ (L 2 ), then u n £ C(R+;L 2 ), by view of Energy Identity ([TIT]) , 
Therefore, it reduces to get a priori estimates of V7r", /i n_1 in L\ oc (£?* r ) (L 2 ), by use of Ct and U n . 

In fact, it will immediately follow by observing the estimates in the demonstration below of the inductive 
estimates for u n and V7r n . 

In the following, we a priori demonstrate inductive estimates for u n and V7r™. The idea is that, by 
terms of r and 

we are going to bound the following quantities in the following order: 

- IIV^IUx.^) - ||V7r"|| z ^ (B .-i) - UAttI^.-!), 



i i . 



which include informations on both low frequencies ||V7r™|| i ^( iP ) and high frequencies ||A7r n ||ji (s s_1 ) 
One also notices that if p > 2, then |V7r™|| i; i^ ^ L 2j readily offers bound on ||A_iV7r™|| i i^ ( - i p' ) . 

(i) First of all, by view of Equation (| 1 1 5[) , , Lemma 14.11 ensures that 

ll« B |lz~ (B . r ) < Ce cwn -^ (\\u a \\ B s r + lift"" 1 - A*V7r1 Zi(B . ir) ) , 

where 

W 1 - 1 ^) := / \\u n - 1 -K n - 1 {p n )- 1 Vp n \\ B * ■ 



o 



By virtue of inductive assumptions (jTTj) for 1 , g n and (f72"j) for it" 1 , we easily derive 

W- 1 ^*), \\h n - l \\~ L ^ {Bir) <CC M r. 
Thus (as ||A™ — 1 II Z°° (_b s ) — ^ Cm) for sufficiently small parameter r, we have 

(127) IK%~( B|r ) <c(\\u4 B s >t + CCmt + CCm\\V^z^ {bi J) < CCj* (l + IP 
(ii) Secondly, inductive assumptions (1721) and (|119[) imply 

ll/i"- 1 !!^^) <cc Ej b t. 

Hence, by view of || Vu" || L ^ {Lao} < {T*)^ 2 \\u n \\ L ^ {B} r) and Estimate (|T27j) . if (T*) 1 / 2 < r, then 
Equation (|12ip entails 

(128) IIVttI^^ < CC E E {T + T\\u n \\ L ~ {B > r) ) < CC E C M E r(l + U n y 
Correspondingly, ||d t u™|| £ i^ ( L 2) is bounded also by above, with some change of the constant C. 
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(iii) Next we also want to bound ||V7T n [|£i ^ (£j>) for p e (1,2), which controls the low frequency. It relies 
on Equation (|122p . Firstly, since 

A&tp"- 1 ), V 2 a(p B - 1 )6L^( B «-i)^4,( L P*) ) with p *^J^> 2 , 

F ' 2 — p 

we have furthermore 

W^Wl^l,) <CC e E t. 

Similarly, we have 

\\(p n - l)d t u n \\ L ^ {LP) , \\{p n u n ~ l - K "- x Vp") • Vu n \\ LULP) < CC E C 2 M E Q r(l + IP). 
Hence Equation (|122p implies, for p £ (1, 2), 

(129) l|V7r"|| L ^ (LP) < CC e C 2 m E q t(1 + IP) < C n (l + IP)r, 

with notation Cn denoting some constant depending on s, d,p, p*, p* ,Ce,Cm , Eq, to be precisely 
determined later. 

(iv) Now, one observes that Estimate (|129p for V7r" implies moreover 

< ^liv^llit^)!!^"!!^,.) ^ + n> 1/s , p e (1,2), 

with some change of constant Cn- On the other hand, if p > 2, then Embedding L 2 (M. d ) «-4 
<^ d/2 (R d ) and (555) also ensure similar interpolation inequality, for any 77 6 [0, 1): 



<C||V7r"|| L ? (L2) (IP) 
< Cna + IP)^. 

(v) At last, Equation (|123p and Estimates (|124p and (I125[) ensure for some 7/ S (0, 1) (still with some 
appropriated constant Cn) 

ll^llz^B- 1 ) < C °m II Vtt" || £ ^ (B| - 1+ , ) +CC M r + CC M r 2 < C n (l + IP)t^, 
which, together with (|120p . (|128p . (|129p and the definition of II™, implies, for r and T* small enough, 

\\^ n h, T , {Bir) <T^. 

Therefore by virtue of Estimates (|127p and (|128p , inductive assumption (j?2")) and (|118p for u™ and 
V7r™ follow respectively. 

(vi) From above, Energy Identity (|117p holds and hence we have 

< C(\\u \\ L , + H/i"- 1 ^!,^)) < \c E E Q . 
Inductive assumption (|119l) is then verified. 
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(130) 



6.2.2 Convergence Part 

Let us turn to establish that the above sequence converges to the solution. As in Subsection 14.2.21 we 
introduce the difference sequence 

{6g n , Su n : VJtt") = (g n - g n ~\ u 11 - u n ~\ Vtt" - Vtt"- 1 ), n > 1. 

When n > 2, it verifies the following system: 

d t Sg n + u n ~ 1 -V5g n - div (k^ 1 VSg n ) = F"" 1 , 
d t Su n + (m"- 1 - K™" 1 Vlogp") • VSu n + X n VSir n = HI 1 - 1 , 
div(5u" = 0, 
(5 Q n ,5u n )\t=o = (0,0), 

where 

F 71 - 1 = -8u n - 1 ■ Vg n - 1 + div (Sk^V g n ' v ), 

if"- 1 = Sh"- 1 - [5u n ~ l - 8n n - l \7 log p n - n n - 2 \75{\ogp) n ) • Vu"" 1 - 5\ n Vir n -\ 

with 

Sk 11 - 1 = n n ~ l - K n ~ 2 , SK 1 - 1 = h"- 1 - h n ~ 2 , 5{\ogp) n = logp n -logp™" 1 , SX n = A" - A n_1 . 

Firstly, since in the case p S (1, 4] we have the embedding B^ r ^ B 1 i \ d ^ i , we just have to establish 

that {(g n , u n , V7r n )}„>o is a Cauchy sequence in the functional space E^(T*) (see (flT)) for the def- 
inition). In fact, we just have to follow exactly Subsection 14.2.21 with some necessary changes in the 
coefficients and source terms switching from System (|78l) to System (|130|) . 

If p is big enough and we consider the limit case (s,p, r) = (1, oo, 1), then it's no more true that 
{(g n , u n , V7r n )}„>o is a Cauchy sequence in i(T*). Indeed, applying "div" to Equation (1130[) .-, entails 

(131) div(A"Vfe n ) = dWH^ 1 - div((u™- 1 - n™- 1 \7 log p n ) ■ \75u n ). 

As presented in the proof of Proposition 14.61 in Subsection 14.1.21 (see (l68|) for details), the following 
estimate 

HdivK" 1 • V*u")|| H */„-i < II^^L./p+iUV^IL./p-i 

p, 1 p , 1 p , 1 

only holds for p < oo. Thus we cannot have V5ir n £ (B^, x ) in general. 

Therefore in the general case, we have to consider the difference sequence in the energy space. One 
wants H™- 1 in L\,,{L 2 ). One pays attention to the following terms in 

(p^^ASb^Va 71 - 1 and (p™" 1 )" 1 Vb"" 2 • V 2 <Ja n_1 . 

We only have Va" _1 , Vfe" -1 in (L°°), and thus one requires A<56 n_1 , V 2 (5a"~ 1 in L\, (L 2 ) and hence 
<5e n_1 in L^,(H 2 ). 

Firstly, it is easy to see that F 11-1 G L 2 oc (L 2 ) and hence taking L 2 inner product between Equation 
(|130p i and 5g n gives controls on 8g n by 8g n - 1 , (5u™ _1 , with small coefficient if restricted on the small 
time interval [0,T*]. 

Next, thanks to the null initial datum for 8g n , by taking derivation of Equation (|130[) . we expect to get 
energy estimates for VSg n . In fact, the equation for VSg n , n > 2 reads 

(132) d t V5g n + u^ 1 ■ V 2 8g n - div (/^V 2 ^") = -V5g n ■ Vu™" 1 + div {V5g n ® Vk"" 1 ) + VF"" 1 . 

The first two terms of the right-hand side are of lower order, while the third one is in L 2 oc (7J _1 ), thus 
taking L 2 inner product works. Therefore, Sg n £ L'^ C {H 1 ) nL 2 oc (H 2 ) ensures Sg n s C(R + ;L 2 ) and 
H™' 1 £ Lf oc (£ 2 ). Thus energy inequality for 5it ra also follows and its energy is bounded in terms of 
Sg n ~ 1 ,Sg n ,6u n - 1 . 

Thanks to the small time T*, we thus can demonstrate that {(p n ,u n )} n >i is a Cauchy sequence in 
C([0,T*];i 2 ) and hence converges to some unique limit (p,u). Furthermore \\p — p n \\ L ^ m^nL 2 » (h 2 ) 
goes to as n goes to oo. Then by use of the high regularity of the solution sequence, we can show that 
(p, u) is a solution. 



44 



Inviscid Low-Mach Number Limit System 



Now we begin to make the above analysis in detail. 
Our goal is to demonstrate that {p n — po} n , {u n } n are Cauchy sequences in C([0,T*];L 2 ) and the limit 
really solves System (|114p and System (|115p . 

Since 5g n £ E (see (I126[) for definition), we can take L 2 (R d ) inner product between Equation (|130[) 1 
and 5g n , n > 2, entailing 

ld_ 

2dt 



\Sg n \ 2 + K n - 1 |V5e n r = - / <fa • V/- 1 ^" - / SK n - L Vp n - 1 ■ V5q 



Thus integration in time and Young's Inequality give 

(133) < CTdl^-^U- (La) + WSu^ 1 ^^). 

By view of the analysis above, Sg n s C(M. + ;L 2 ) n L 2 oc (H 2 ), n > 2, (notice that it is not clear that 
Sg 1 = p 1 - p € £j oc (H 2 )). So we can still take L 2 (R d ) inner product between Equation (|132[) and V(5g": 



l_d 
2 eft 



f \VSg n \ 2 + [ K n - 1 \V 2 5g n \ 2 = - [vSg n -Vu n - 1 -VSg n + VSQ n -V 2 Sg n -VK n - x + F n - 1 AdQ n . 

JR d JM d J 



Integrating in time also implies 



< 



\\ySg n \\ L ^ {L2) + \\V 2 5g n \\ L ^ {L2) 

c(||v« n - 1 || i? , (i o. ) ||v^ n |U» (L2) + ||Vp"- 1 || i; , (i o. ) ||v^ n |U ? o + IIF-- 1 !!^^)). 

Since 

II^-Il^) ^^II^^II^.^+tII^-I^^j + CmIIv^"- 1 !!^^)), 

for the above small r and T*, 

[|W|| £ ~ (£a) + ||v 2 ^"||^, ( l=) < CriKJe"- 1 ,^"- 1 )^- (La) +ccw[|V5e n - 1 [[ i3 ,. (£a) . 

We can substitute (| 133[) into the term || VSg"^ 1 \\ L 2 ( ^2) above, and then sum up these two inequalities, 
entailing 

(134) \\5 Q n \\ L ~ + IIWL^hx) < CCmtIK^"- 1 ,^"- 2 ,^"- 1 ,^"- 2 )!!^, (L2) . 

Now we turn to <5w n . We rewrite 5h n ^ 1 as 
— ^-(A^^Va™- 1 + A6"- 2 V(5a n - 1 + cV 1 - 1 • VV" 1 + u n - 2 • V^a"" 1 + VStf 1 ' 1 ■ VV" 1 

pn-l V 

+W l ~ 2 • V^a"" 1 ) + ((p" -1 ) -1 - (p'^r 1 ) (Ab n - 2 ■ Va n - 2 + u n - 2 ■ V 2 a n - 2 + V&"~ 2 • V 2 a™~ 2 ). 

By virtue of WASb 71 ' 1 ]]^ < C||5e n_1 [|H= and 

HA^- 1 !!^^) < C|||V P "- 1 | 2 || L|t(L=0) + CHAp"- 1 !!^.^) < CC m t, 

we have also from the above inductive estimates that 

WSh"- 1 ^^ < CCuTiWSg"- 1 ^^ + llfo"" 1 !!^^)) +CC E EoT\\5g n - 1 \\ L ~ (L2) . 
Similarly, 

WHr'h^iL*) < C(C M + C E £ )r(||^ n - 1 |U ?r ,(H=) + ||5n"- 1 |U ?f ,( L2) )+Cr 1 / 2 ^||^"|| i? o (i2) . 
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By view of Equation (|114p and div Su n — 0, taking the L 2 (R d ) inner product between Equation (I130IU 
and p n u n and then integrating on the time interval [0,T*] yelds 

(135) \\5u n \\ L ~ iL2) < CHlfr 1 ll^(^)- 

Combining Estimate (|134l) and (1135[) entails, for sufficiently small t, depending only on the space 
dimension d 7 on the indices (s,p, r) and on the constants Cm ,Ce, Eq, 

\\5e n \\L~(m)nL* T ,( H z) + \\5u n L? t (L*) < |-|| (.5^- 1 , 2 , ,5^- 3 , 1 , 2 , ) | c ^ 2) . 

Thus \\(Sg n , 8u n )\\ L cv t ( L 2^ converges. Since Sg n £ C(R + ;L 2 ), the Cauchy sequences {g n } and {u n } 
converge to g and u in C([0, T*}; L 2 ) respectively. It is also easy to see that 

Y \\ s e n \\L~(Hi)nL 2 T ,(H2), Y ll^ll^^i 2 )' Y ll^"^ 1 !!^^ 2 ) < +°°- 

n>2 n>2 n>2 

Writing 

div ((/-' - k"- 1 Vlogp") • VSu n ) = div (Su n ■ V(u n_1 - re™- 1 Vlogp") + <5u"div (k^V log />")), 
from Equation (|13ip we get 

l|vfei^ (i2) < c(||iC _1 IUi,.(£») + CM||«s« n lU» (£a) ). 

Thus 5^2° ^^"""lli^,, (L 2 ) a l so converges and hence V7r ra converges to the unique limit V7r in L^.(L 2 ). 

Now one notices that, by interpolation between L 2 and i?p r , «" converges to u in (7([0, T*]; j 2 ), 
for instance, and — K T * _:l (/j n ) — 1 Vp n is at least in Lf oc (B\^). It entails 

(u n - 1 -/s n - 1 (p n )- 1 Vp n )-Vu n -(u-«pVp)-Vu-»0 in l£.(.B^ 4 ). 
Thus a direct computation ensures that (g, u, Vir) solves System ([7]) and is in E* r (T*) by Fatou property. 

6.2.3 Uniqueness Part 

The proof of uniqueness just follows the idea of the convergence part. More precisely, as in Subsection 
14.31 take two solutions (pi, u\, V7Ti), (p 2) tt2> Vt^) £ E£ r (T*) of System ([7]) with the same initial data, 
such that pi — l,p2 — ljiti, u 2 G L^ t (L 2 ), Vpi,V/02 £ £y,(L 2 ). Then the difference (<5p, (5m, V<57r) = 
(pi — P2> Wi — w 2, V7Ti — V7r 2 ) verifies 

+ ui • V5g — div (KiV<5g) = —5u ■ \7g 2 + div — K 2 )Vg 2 ), 
<9 f <5u + (tti + V6i) • VJu + AiV^vr = hi-h 2 - {Su + VSb) ■ \7u 2 - 6AV7r 2 , 

div 5u = 0, 
(5ft5u)|t=o = (0,0), 

with the notation ft, = and analogous for hi, Xi and ftj. 

Similarly as Convergence Part, we can get 

¥p\\l?{li) + IIVMI^) < C(\\Vp 2 \\ L 2 {L ^ ) \\Sp\\ Lr(L 2 ) + \\Vp 2 \\ L i {L ~ ) \\6u\\ Lr{L 2 ) ), 

||VJp|| L -(L») + l|V 2 <M| i?(i2) < c((l|Vui|| L i (L oo) + llVKilliKioojJHVJplli-.^) 

+ (II A P2||l 2 (L~) + l|Vpi|| i 2 (i0 o ) ||Vp2||Lf>(i«>))||5p|U«»(ia) 

+ IIVp 2 ||L 2 (L~)IIHUr(£ 2 ) + II v piIU?°(£~)II V(5 pIIl 2 (l2)), 

and some similar estimate for ||^u||l°°(l2), which we omit here. Thus, on sufficiently small interval [0,i], 
dp = 5u = 0, the uniqueness holds. Then we recover uniqueness on the whole existence time interval 
[0, T*} by use of classical arguments. 
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6.3 Remark on the lifespan in the 2-D case 

In this subsection, we want to give a better estimate for the lower bound for the lifespan of the solution 
in the case of dimension d = 2. 

As mentioned before, the global-in-time existence issue in dimension d = 2 for the classical homoge- 
neous Euler system (i.e. p = ~p constant in system 0) has been well-known (see [?] by Wolibner, or also 
Chapter 7 of [?] and the references therein for a survey on this topic) . The key to the proof is the fact 
that, if we define the vorticity of the fluid as 

(136) u := diu 2 - &2U 1 , 

then this quantity is conserved along the trajectories of the fluid particles, i.e. it fulfills the free transport 
equation 

(V) d t uj + u ■ Vuj = . 

For non-homogeneous perfect fluids, see system (fTTj) . the previous relation (V) is no more true, due 
to a density term which comes into play combined with the pressure. Hence, it's not clear if solutions to 
(fTTj) exist globally in time. However, in [?] it's proved that, for initial densities close to a constant state, 
the lifespan of the corresponding solutions tends to +oo. 

Theorem 12.111 gives us the analogous result for our model. The idea is to resort to the vorticity in 
order to control the high frequencies of the velocity field. The vorticity uj of the fluid is still defined by 
formula (|136p . where u solves ([7]). However, it's easy to see that actually 

uj = d\v 2 — ; 

now, from (0]) 2 the equation for uj immediately follows: 

(137) d t uj + v ■ Vcj + uj Ab + VA A Vn = , 

where we have set VA A Vn = diAc^n — c^Ac^n. 

In order to bound the vorticity, it's fundamental to take advantage of a new version of refined estimates 
for transport equations in borderline Besov spaces of the type x , proved first by Vishik in [?] and then 
generalized by Hmidi and Keraani in [?] (see also Chapter 3 of [?]). They state that the B^ x norm of the 
vorticity grows linearly (and not exponentially, as in the general case) with respect to the Lipschitz norm 
of the solenoidal velocity field. Even if here we don't have the divergence-free condition for the transport 
velocity v, the proof in [?] still works, except that one asks for an additional regularity on divi>. More 
precisely, the following lemma holds. 



Lemma 6.6. Let us consider the following linear transport equation: 
(138) 



dtuj + v ■ \7uj = g, 
uj\ t =o = uj . 

For any /3 > 0, there exists a constant C depending only on d, ft such that the following a priori estimate 
holds true: 

(139) 11^)11^ < c([[w || B D >ii + y| L|(i3 c o i) )(i + v(i)), 

with 

V(t) := / ||V»||i- + ||divw|| B , dt 1 . 



Proof. We will follow the proof of [?]. Firslty we can write the solution uj of the transport equation (|138[) 
as a sum: uj — X)fc>-i ^k, with uj^ satisfying 



(140) 



d t uj k + v ■ Vuj k = A k g, 
uJk\t=o — AkUJo- 
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We obviously have from above that 

(141) lhk(*)lk~ < HAjfcWoHi-. + / \\A k g\\ L ~ dt'. 

Jo 

By Proposition 13. 121 for any e G (0, 1), we have the following a priori estimate in Besov space x : 

(142) WMQWb^ < (\\A k cj \\ Blo i + \\A k g\\ Li(BL i) ) exp(c|| Vu^i^-)) . 

In order to get a priori estimate in Besov space B~^ 1 , after applying the operator A.,- to Equation 
(| 140(1 . we write the commutator [v, Aj] ■ Vui k as follows (recalling Bony's decomposition ([50)) and denoting 

v := v — A_it>) 

[7~, AjYVLOk+TA^^v+RiAjVLOk^-AjlJv^-Ajdiv {R(u k , v))+A j R(u) k , div v) + [A_it>, A^-Vwfc. 

Then, V/3 > e, the L 00 -norm of all the above terms can be bounded by (for some nonnegative sequence 
ll(ci)ll* = l)= 

C(d,p)2-^c j V(t)\\u k \\ B - :y 
Thus, we have the following a priori estimate in the space B^^. 

(143) HwfcWIlB-., < (||A^o|| B - :i + \\A k g\\ Ll{B -*^ exp(cV(*)). 

On the other side, one has the following, for some positive integer N to be determined hereafter: 

j,fc>— 1 |i-fc|<AT | j — fc | > JV 

Estimate f| 141[) implies 

£ IIA^wjblli- < nJ2(\\ A ^o\\l^ + \\A k g\\ Ll(L ^ < iyr(||w ||flo iil + IMUk^))' 

|j-fc|<AT fe 

while Estimates (|142l) and (I143[) entail the following (for some nonnegative sequence (cy) e f 1 ): 

HA^fclli- <2- £ l fc -^(||A fcWo || L oo + ||A fe5 || L j (i oo ) )exp((7V(i)), 
which issues immediately 

]T IIA^kHi- < 2-^([[w || B o )ii + hhliB^)) exp(cV(t)). 

b-fe|>iV 

Therefore, for any /3 > 0, we can choose e G (0, 1) and JVeN such that e < f3 and TVe log 2 — 1 + CV(£). 
Thus the lemma follows from the above estimates. □ 

With the above lemma in hand, we can now begin to prove Theorem l2.11l In view of the continuation 
criterion, without loss of generality, we will always assume in the sequel 

s = 1, p = oo, r = 1. 

In what follows, we will resort to the notation of Subsection 15.21 That is, we set 

R(t) = llellifCB^), s ( t ) = II^N^L,!)' = NUr^i,,!)' 

and the time as defined by (|57|l . 
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First of all, let us consider the density term. Under our hypothesis, we get the energy equality (|24[) . 
Moreover, arguing as in the previous subsection, we recover estimate ([88)1 in the time interval [0,Tr]. 

Now, let us consider the velocity field. First of all, let us summerize the following inequalities for the 
nonlinear terms, which will be frequently used in the sequel: 

(144) ^Hp^b^ZMb^ZMb^^S'-, 

(145) ||A6Va|| L2 < p]]^ ||Va|| i2 < \\g\\ Blo jVp\\ L 2 < R 1 / 2 ^') 1 ' 2 ^ p\\ L 2 < R\\Vpf L2 + S'; 

(146) ||(u + V6) • Vu\\ L 2 < ||V«||i-(||Vp[| ia + \\u\\ L2 ) < U(\\V P \\ L 2 + \\u\\ L 2). 

Similarly as the above inequality (|145p . one has 

(147) ||V&-V 2 a|| L 2 <R\\Vp\\ 2 L3 +S', \\u ■ V 2 a\\ L 2 < R\\u\\ 2 L2 + S' . 

It is time to bound the velocity field u by use of the above inequalities. Firstly, by separating low and 
high frequencies of u and by use of Bernstein's inequalties and of a Fourier multiplier of order —1, we get 

(148) U(t) < C (\\u\\ L * + IMIflO^) . 

In order to handle the energy of the velocity field, one takes the L 2 scalar product between (TTUJ) and 
u and performing standard computations (recall also ([25)1 and the following arguments), getting 

K*)||i» < c(\\u \\ L 2 + J ||div (v®Va)\\ L 2dT 
Therefore, due to Inequalities (| 145)) and (|147[) . it follows that 

(149) \\u(t)\\ L 2 < C (\\u \\ L2 + J*(R(\\Vp\\ 2 L 2 + \\u\\ 2 L 2) + S') dr^j . 
Now, applying Lemma 16.61 with f3 = 1 to Equation (|137p , we find 

\\uj{t)\\ BLi <c(\\uj \\ Blai + ^ f |vAAVn + W A&| so dr) (l + J (|| Vu\\ L ~ + \\ V 2 b\\ BL J drY 

By use of Bony's paraproduct decomposition (see also Section 4.2 of [?] for the first inequality) 

[[VAAVttUb^ < HVplls^ [[Vllllflo^ 
|[o;A6|| B o i < \\u\\ B o \\Ab\\ Bloi . 

OO , 1 00,1 OO , 1 

Hence, by virtue of ||w|| B o < U and Inequality ()144p . we immediately gather 

(150) HwCOlU^ x <c(u Q + J (iJHvnn^ x + us'} drj (1 + ||v«|| £ -odT + s\ . 

It remains to deal with the pressure term. First of all, we have 

(151) WWWb^ < IIVttIIbo^ + II^VoHb^ a < HVvrll^^ + HqW^. 
Thanks to equation ([7]) , , and by use of Proposition 13.101 we get 

\\dtg\\ BLl < \\u-v Q \\ BLi + \\kVq\\ b . Li < US' + S'. 

Let us now focus on ||V7r|| B o ^. Actually, we will bound the l norm, as it's not clear for us how to 
get advantage of the weaker norm in (|150|) . The analysis is mostly the same performed in the previous 



49 



F. FANELLI and X. LIAO 



subsection for the general case: so we cut low and high frequency, and we are reconducted to consider 
||V7r||i,2 and ||A7r|| B o . Finally, as 6 > 1, we have 

||V7r|| Bj)ii < C {l + R s ) (HVttII^ + ||div(«.V«)|| B2eii + \\h\\ BLi ), 

with the controls (by what we established in the case of higher dimension) 

\\HbI o1 £ (1 + R)US' + (l + R)RS' 

and, as div (u ■ Vu) = J2i j ^Tg jU idiU 3 + djR{u l ,8^) (thanks to the divergence-free condition over u) 
and analogous for div (V6 • V«), 

|[div(«- V«)|| B o )i < U 2 + US' . 

On the other hand, observing that ||/i||l2 < ||div (v Va) H^a and using Inequalities (|145l) . (|146|) and 
(|147p , one has also 

|]V7r|| i2 < ||A|| ia + || V • Vu|| L2 < R{\\Vp\\ 2 L2 + \\u\\l 2 ) + S> + U(\\V P \\ L , + \\u\\ L2 ). 
Let us define 

X(t) := U(t) + \\u(t)\\ L 2 = \\u(t)\\ L2nBloA . 

So we get 

(152) ||Vtt|| b i < C (l + R 6+2 ) f||Vp||| 2 + S' + X 2 + XS' 



and, by (|15ip , the same holds true also for [|VII|| B o . 

Therefore, Estimate (I150[) for the vorticity becomes (denoting Xq = X(0)) 

IIwWHbo, 1 < C ( X + S + Xdr^j (^X Q + J q ( 1 + rS+3 ) II Vp||i a + R S' + R X 2 + X S'^j dr^j . 

It's now time to insert the above estimate and (|149[) into (|148|) . Keeping in mind that, in [0, Tr], 
(153) R(t) + S(t) < CR cxp (° ( 1 + x ( T ))d^j , 

we finally find 
(154) 

X{t) <c(l + S + j* Xdr^j lx + Ro(l + R s + 3 )e c ^ 1+x ^ J* ||Vp||| 2 + 



+ i?o(l + < f4 )e C/ ° (1+A ' ) +R a (l + R 5 +3 )e c fo(^) / X 2 + (1 + R s + 3 )e c K ^ / XS' dr 



r 3 r 4 
We define Tx as the following quantity (with Ti, i = 1, • • • ,4 defined as above) 

T x := sup{i | Tx(t) < \\qo\\ 2 L 2, r a (i) < 1, T 3 (t), T 4 (t) < 1 + \\g\\ 2 L2 + X }. 

Then, noticing S < T2, one easily arrives at the following bound for X(t) with t in [0, Tr] n [0, Tx] (with 
some positive constant still denoted by C) 

X(t) < C(l + ||fio||^+X ) (l+ f X(r)dr 
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Hence, set Tq := C(l + H0OHI2 + Xo), then by Gronwall's lemma we get 
(155) X(t) < T e Tot , 

and the norm of the solution can be controlled by the norm of the initial data only. 

Our next task, in order to complete the argument, is then to prove that T, defined by (|29[) with 
yjl = 5 + 4 and small enough constant L, is smaller than both Tr and Tx- First of all, thanks to (|155D , 
for t e [0,T X ] n [0,Tr], we have 

fx<e T *\ fl 2 <V» f , f\l + X)<2e r ° t , e c & 1+x > < e 2C <?°\ R + S < CRae 20 * * . 
Jo Jo 2 Jo 

With the above bounds in hand, one just has to show 

r x (T) < \\Qof L2 , T 2 (T)<1, T 3 (T), r 4 (T) < l + ||e ||| 2 +X and f i? 3 <2i? . 

Jo 

We will check the above bounds one by one. 
By use of the energy equality (|24[) for the density, it's then easy to see that Ti(T) < \\qo\\ 2 L 2, with T 
defined by (f2l?)) . It is also easy to find that ^(T) < 1. Now noticing that (<r < e CT ) 

r 3 (t) < Ro(l + R 5 +3 )e 2Cerot §e 2rot < ^i?o(l + ^ +3 )e 3Ce2rot , 
we hence have ^^(T) < 1 + ||£o||f,2 + Xq. Similarly, since one has 

r 4 (t) < (l + i?^ +3 )e 2Cerot r oe rot /' s' < r (i + R^ 3 )e 5CeTot CR , 

Jo 

thus T^(t) < 1 + || go II 2 2 + Xq. Finally, it's not hard to check that 

C R 3 <C 3 R 3 e 6Cer ° T <2R . 
Jo 

This completes the proof of the theorem. 



6.4 The case when the density is near 1 

Let us give here a sketch of the proof to Theorem 12. 121 In fact, we will establish only a priori estimates, 
the rest of the proof to existence and uniqueness being similar to that of Theorem 12.11 

So, recall that p G (1, +oo) and ||go||s= r — c ' wnere c is a small positive constant. Moreover we will 
use the convention that r = lifs = l + d/p. 

We first focus on the density equation §7$ 1 . We start from inequality (|5^|) . in which we use also the 
second relation of (|54[) to bound the first commutator term. 

We apply instead (|4"o| to control the second commutator term: with a = s, 6 = rj = 1/2 and e = 1/2, 
it entails 



2-" 



\\Lp 



dr < C 



livell 



I ells;, dr 



1 



llellz 



Putting all these relations together, we get 
for constant C depending only on indices d, s,p, r and on p„, p* . 



iVul 



•C||Ve|||- dr 
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Now, if we take T €]0,Tb] small enough, such that, for instance, 
(156) exp^C< + C^ ||Vu|| B .-i +C\\Vg\\ 2 L oodT\ <log2 : 



then the previous inequality becomes 

(157) llellz-fB. 



I^HiJ(B|+ 2 ) < 2C*c. 



Let us now consider Equations (JTJj and ([7]) 3 for the velocity field and the pressure term. Going along 
the lines of the proof to Proposition 14.61 we can see that ((62)) and (j63|) still hold true, with the transport 
velocity w being u + Vb. 

Here we handle Vir as done in the previous section, that is, we make use of (|120[) . Firstly, notice that 
the equation for n can be rewritten as 



(158) 

which (formally) implies 



An = div | h - (u + Vb) ■ Vu + -Vtt j . 



Vir = V (-A) -1 div I - h + (it + Vb) ■ Vu - - Vtt 

P 



Hence, by Calderon-Zygmund theory and (I120|) we infer (with w = U + Vb) 
(159) ||V7r|| Zi(fl . j < C [\\h\\ ll(BSp r) + Wd^d^W^-^ + \\w ■ Vu\\ Ltm 



e -v* 



The term \\diw3 djii l \\ z^b^ 1 ) can controlled as done in For \\w ■ Vu\\ L i^ LP ^ we notice that 



\\w-vu\\ L]{LV) < / |HU*IIVtt|U- < / (||v«|U- + \\Vb\\ LP )\\u\\ B}r 

Jo Jo 



Finally, we have 



Q - Vtt 



< c Mz Ti B^ \Mz U B l:r) 



Therefore, putting all these inequalities together into (|159p and using (|157[) . we finally arrive, for c 
small enough, at the following relation: for all t £ [0, T], 



(160) 



l|V7r|| Z x (B ^) < Cl\\h\\ Zi(B . r) + / (||V U || B .-i + ||Vp|| B . ir )|| U || B . r dr 



Now, the analysis of the nonlinear term h can be performed as before, and, up to taking a smaller T 
than the one defined in (|156|) . this allows us to close the estimates on some suitable time interval [0,T], 
for sufficiently small c. 



A Appendix 

First of all, let us recall an easy version of Young inequality: for all 9 <e]0, 1[, all e > and all real 
numbers a, b > 0, we have 

(1) ab < 6 e^ 1 -^ a 1 ' 8 + (1 - 9) e fcVd-fl) . 

It is the main ingredient to the proof of Lemmas 14.21 and 14.31 
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Proof of Lemma \4-S\ We decompose the commutator by use of Bony's paraproduct: 

(2) [^,A 3 W = R)M) + RjM) + + + R 5 jM) 

where, setting yp — (Id — A_i)<^, we have defined 



R)M) 
R 2 iM) 



T kjV^ = z2 s k+2A j Vtj} A k ip 
k 

—AjTy^ 

[A_i^,Aj] VV>. 



Let us point out here that, in the following, in general we will derive two types of estimates for each term, 
corresponding to (|44[) and (|46[) respectively. For simplicity, in proving (|46p we will consider only the case 
6 = rj\ the general result follows from easy changes in the proof. 

In the following, for the notational simplicity, d x will always refer to one of the partial derivatives d xi , 
/ = !,•■• , d. Let us first consider the term 



K 1 



2" 



\d x R]\\ LP dr 



where one finds already 



h{z) z ■ VSv-wix - 2~ J \z) A u Vtf>(x - 2~ 3 z) dX dz 



As it is supported in a ball of radius 2 3 , the derivative d x gives (by Bernstein's inequalities) a factor 2 3 . 
Hence we get 



n 1 < 



1/r 



we can decompose the terms in the integral in the following way: 

V" llV^-i^ll^ ||A,V^|| LP < 2^ HVSf-x^ll^ HA^IlL 7?W* ||A^||£ e . 

Now, we apply Young inequality |1| to separate the two factors; therefore, using Minkowski inequality 
for the first term, we get, for some constant C: 



(3) 



TZ 1 < 



a 



r(i-e)/e 



I|V^||^||V||b-^+ = II 



Let us now handle 









K 2 := 




< 




Jo 





2 iCT / V VA^IU- ||A^|| iP 

J ° /»>J-2 

2"/ £ 11^+2 VA^Hioo ||^A^|| L pdr 



One notices that, since /i > j — 2, we have 

VA^||i=o ||A^|| iP < 2"||VA i ^|U»||A^|| iJ , < ||VA^||i=o HVA^HiP 
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Hence 



n 2 < 



< 



V° [ Yl ||VA^|k~||VA^||^dr 
[ ||V^|U- E 2^^2^||VA^|| iP dr 



On one side, we just do exactly as above (the way to obtain (O ): if a > 0, then we have 

C0 



(4) 



(i-e)/e 



On the other side, by the relationship 2 3 * ||Aj-/||x,oo < ||Aj/||lp and Young's inequality, we further 
derive 



n 2 < 



< 



< 



/■* 2 i(-+i-f)|| Ai¥) | |ioo £ 2 0'-M)| ( 2 ^||VA^|| iP ) dr 

( ct+1 )||A^||l p 2^||VA^|| LP 
£ 2^)fc,||V^|L, - 2i(^+( 1 - e )-)||A,^||i r(iJ)) ||A^|| i 7 ( % 



(5) 



C6» 



{cj). g r 



i7(Bp P oo) 



Moreover, since 











ft 3 := 




< 


2J( CT +!) /" 




Jo 




Jo 



2^+i) f J2 ws^MUoo \\A^\\ LP 



we can immediately see that ^ holds also for 1Z 3 . In order to get the form like JS}, we consider the two 
cases a < | and a > ^ separately. Here we notice that if a = ~ and r = 1, then taking the integral in 
time and the £ r norm can commutate and hence the classical result for commutators (see Chapter 2 of 
[?]) gives this lemma. Now ifc^ < 1 + -, i = 1,2, then 



K 3 < 



■e)Oa-l-#)| 



which gives (JSJ for 1Z 3 . On the other hand, &{ > 1 + ~, i = 1, 2, then as above 



ft 3 < 



o 

C6> 



HWIU- (2^ +1 )||A^|| 



(G) 



^■^=m^U fRr+1 JIWIIl.- ( l-) + (i-flj^iiv^o.) 

ij t (-Dp,,- ) 
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where the last inequality follows by the embeddings Lj(L°°) L\{B p l) Lj(Bp 2 r x ). 
The same holds true also for 



TZ 4 := 



2JCT / \\ d * R j\\ L * dT 



Let us first sketch the case a > -. We can bound TZ 4 in the following way: 



n 4 < 



M>i-2 



Hence (jH and © follow immediately for TZ 4 . If a < ^, then we consider two cases p < 2 and p > 2 
separately, which both give ©. In fact, if p < 2, then p' > p and we have for cr + 1 > —jr, 



TZ 4 < 



< 



2^ +1+ ^ ) ||i?4|| Ll 



/j>j-3 



- 1+ F)(2-(f+ 1 )||A^|| iP ) (2 
which gives ([5]). Otherwise if p > 2, then we have for cr + 1 > — ~, 

JO f 

/* £ 2^( CT+1 +f)(2^f +1 )||A^|| LP )(2^||AMV^||Lp) 



ft 4 < 



< 



H>j-3 



which yields ([3]) also. 
Finally, last term 



TV := 



2" / II^^H^ dr 



can be handled as TZ 1 , leading us to the same estimate as ([3]) and so to the end of the proof. 
Proof of Lemma \J73[ We have to estimate 

1/r 

/ f r t \ T ) 



□ 



^ J\\^{f g)\\ LP dr 



Using Bony's paraproduct decomposition, we can write 

(7) A..:/,,; = A, 7}// + A//,/ + A, /i' ;/.//:■ . 

Let us consider the first term in the right-hand side of the previous relation. Due to spectral local- 
ization, we infer that there exists a positive integer lq such that, for all j, 



k-j|<to 



\^T f g\\ LP < C\\f\\ L ~\\A jg \\ LP 
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for a constant C which depends only on to- Now we integrate in time, we apply Young's inequality ([I]) 
and we pass to the l r norm: using also Minkowski's inequality we gather 

H T /5llz; ( B pV ) < ^ J of\\f\\ 1 ^\\9\\ B ^dT + C 2 s(l-e) \\9h i(B;%) ■ 

The term AjT g f can be treated in the same way, and this leads us to an analogous estimate, in which 
however the roles of / and g are reversed. 

Let us now consider the remainder term. Thanks again to spectral localization, we have, for some 
positive integer jo independent of j, 

A 3 R(f,g) = E E A, (A M / A,.g) . 

v>j-Jo \n—v\<l 

Therefore it follows that 

\\R(f,9)\\ Il(B „ r) < C I Y, ( E 2C*-">" 2"" f \\A v f\\L-\\K9\\L»dT 

\j>-l \v>j- M J ° 

Now we apply ([T]) to the product in the integral and Young's inequality for convolutions: we finally arrive 
to 

||i*(/,<?)|| £ x (BpV) < /ll/ll^ \\9\\b^Jt + C 2 e(l-9) \\g\\z UB;2r) , 

and this completes the proof of the lemma. 

Let us point out also that Remark 14.41 easily follows slightly modifying the previous proof. □ 
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